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We investigate the influenee of the thermal properties of the boundaries in turbulent 
Rayleigh-Benard convection on analytical upper bounds on convective heat transport. 
We model imperfectly conducting bounding plates in two ways: using idealized mixed 
thermal boundary conditions of constant Biot number 77, continuously interpolating be- 
tween the previously studied fixed temperature (77 = 0) and fixed flux (77 = 00) cases; 
and by explicitly coupling the evolution equations in the fluid in the Boussinesq approxi- 
mation through temperature and flux continuity to identical upper and lower conducting 
plates. In both cases, we systematically formulate a bounding principle and obtain ex- 
plicit upper bounds on the Nusselt number A'^m in terms of the usual Rayleigh number Ra 
measuring the average temperature drop across the fluid layer, using the "background 
method" developed by Doering and Constantin. In the presence of plates, we find that 
the bounds depend on cr = d/X, where d is the ratio of plate to fluid thickness and A 
is the conductivity ratio, and that the bounding problem may be mapped onto that for 
Biot number 77 = cr. In particular, for each a > 0, for sufficiently large Ra (depending 
on cr) we show that Nu < c{a)R^^^ < CRa^''^, where C is a cr-independent constant, 
and where the control parameter i? is a Rayleigh number defined in terms of the full 
temperature drop across the entire plate-fluid-plate system. In the Ra — > 00 limit, the 
usual fixed temperature assumption is a singular limit of the general bounding problem, 
while fixed flux conditions appear most relevant to the asymptotic Nu-Ra scaling even 
for highly conducting plates. 



1. Introduction 

The Rayleigh-Benard system, in which a fluid layer between two parallel plates is 
heated from below, is a popular model system for the experimental and theoretical in- 
vestigation of the important phenomenon of co nvection, in which dens i ty changes due t o 
heating give rise to buoyan cy- driven fluid flow (|Normand et al\ (|1977^ : iKadanofj (l200lh : 



Cross fc Hohenberd (Il993l) ). With sufficient heating the flow becomes turbulent, and the 



spatiotemporal dynamics become inaccessible to a detailed analytical or experimental 
understanding; instead, one focusses on bulk statistical properties. Of considerable inter- 
est is the dimensionless Nusselt number Nu, which measures the averaged total heat flux 
relative to what it would be in the absence of convection, since the convective fluid mo- 
tion transports heat upward more efficiently than would be achieved by pure conduction 
with the same overall temperature gradient. 

In particular, much research has concentrated on trying to understand the dependence 
of Nu on the (averaged) temperature difference across the plates, represented in nondi- 
mensional form by the Rayleigh number Ra, which measures the relative strength of 
buoyancy and dissipative forces. This dependence is often assumed to take the power 
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law form (with possible logarithmic corrections) Nu ^ CRa^ . Here the prefactor C may 
depend on the geometry of the experimental apparatus (for instance, the aspect ratio of 
a typical cylindrical cell) and/or the Prandtl number Pr. 

In spite of numerous studies over the years, a consensus on the precise form of this 
scaling relationship (especially in the large- i?a limit) remains elusive. Expe riments have 



typica ll y found the exponen t p to lie in the range 1/4-1/3 ( s ee for instance iHeslot et al 



(ll987l): lGla"zier et a/J (( l999'):'Nieme la fc Sreenivas an' ('2006&'): 'Funfschilling et g /J (|2009|). 
or the reviews bv iKadanoff (,200 lE " iProcaccia fc Sreeniva sa n (2008) and .Ahlers et al 



( 20096 )). although higher values have also been reported ( Chavanne et al. ( 200l[ )). Phe- 



nomenological models have also made various predictions, ranging from the early values 
p — 1/3 (.Malkusi (|1954) ) and p — 1/2 (Kraichnan (.1962.)). both suppor t ed by dimen- 
sional arguments, to a more recent model due to iGrossmann fc Lohsd (2000|), which 
predicts different superpositions of scaling exponents i n diffe rent parameter regimes. 
Meanwhile, the numerical investigat ions of Amati et all i 2005 ) in cylindrical geometry 
(performed at higher resolution by IStevens et all ( 2010f) ) found the scaling p = 1/3 
(though, surprisingly , the r ecent two-dimensional, horizontally periodic computations of 
Ijohnston fc Doering (l2009l) are consistent with p — 2/7). 



Effect of imperfectly conducting plates bounding the fluid: 

While some of the variations in the observed results and discrepancies between ex- 
periment and theory may be due to sidewall conductivity, Prandtl number variability, 
non-Boussinesq effects, geometry or other factors, recent attention has increasingly fo- 
cussed on the influence of the thermal properties of the fluid boundaries. The standard 
assumption for Rayleigh-Benard convection is that the upper and lower boundaries of the 
fluid are held at uniform and fixed temperature; this is equivalent to the bounding plates 
being perfectly conducting. In experimental situations, however, the thermal conductivity 
As of the plates is finite, though typically much larger than the conductivity A/ of the 
fluid. In the convective state, the rate at which the fluid transports heat is effectively 
comparable to that which would ensue from conduction with conductivity Nu Xf. Hence, 
for sufficiently strong heating, the assumption that the plates transport heat much more 
efficiently than the fiuid, and are able to maintain the fiuid boundaries at constant tem- 
perature, loses validity; indeed, in the asymptotic high- i?a limit, one might expect that 
relative to the fluid, the plates are effectively insulating. 

A basic consideration in investigating the influence of poorly conducting boundaries 
on convection is the choice of thermal boundary conditions (BCs). Numerous researchers 



have concentrated solely on id ealized fixed flux conditions corresponding to per f ectly insu 
lating boundaries (for instance'Cha pman fc Proctor ( 1980l ): [Otero e t oL Verzicco 



12008); Johnston fc Doering (2009 )), while other studie s (inclu ding[Sparrow et al. (19M); 
[Certsberg fc Sivashinskvl ()l98ll ): IWesterburg fc Bussd (2001)) have imposed more gen- 
eral mixed conditions of flxed Biot number rj at the fluid boundaries. Note, though, that 
the Bio t number in g en eral d epends on the horizontal "disturbance" wave number in the 
plates ( Normand et al 19771 Section V.C.I). For strong driving (high Ra), the temper- 
ature distribution in the plates is unsteady and contains a superposition of horizontal 
wave numbers, so that even mixed, fixed 77 conditions form an approximation to the 
experimentally more realistic situation of a fluid bounded by plates of finite width and 
conductivity. Consequently, some authors have studied the effect of imperfectly conduct- 
ing boundaries by directly incorporating plates in their models, for the study of both 
the con vective insta b ility and the weakly n onlinear behaviour beyond transition (for 
instan ce Imirle et all (Il967t ): IProcto^ (Il98l[ ): Ijenkins fc Proct or' ('1984'): 'Holmedal et al 



fc Sreenivasan 



(|2005l )) and for high-J?a convective turbulence (jChilla et an(i2004 ) : .Verzicca (2004) ). 
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The influence of tlic plate thermal properties on the initial instability of the conductive 
state and the weakly nonlinear dynamics and patteri i formation beyond instability has 



been s t udied intensivel y since t he pioneering works o f ISparrow et ali (|1964 ): iHurle et al 



(|l967t ): iBusse fc RiahU (|l980[ ): IChapman fc Proctoil (jigSOD and others. Their effect on 
heat transport in turbulent convection has, however, only been considered much more 
recently, though it is now receiving attention in the context of experiments, numerical 
computation, phenomenological modelling and rigorous analysis. In the latter category 
is the study by lOtero et al. ( 2002.) , who considered analytical bounds for fixed flux con- 
vection (perfectly insulating boundaries), as discussed further below. 

The suggestion that the finite (even if larg e) heat capacity and conductivity of the 
plates would affect heat transport was made bv lChaumat et al\ (|2002l ). who subsequently 
extended their phenomenological model to propose a criterion for sufficient ideality of 
the plates' thermal properties for the Krai c hnan jo =1/2 "ultimate re g ime" to develop 
(jChilla et al\ (|2004) : see also iRoche et al] (|2005l )): while iHunt et al\ (|2003f) modelled 
the effect of the thermal diffusivity of the lower plate on plume formation and eddy 
motio n. On the basis of extensive numerical studies with varying plate properties lVerzicco 
(|2004l ) concluded that the effects of the plates are governed by the ratio of the thermal 
resistance of the fluid layer to that of the plates, and proposed a model quantifying the 
resultant effect on t he Nusse lt number, w hich was partially conflrmed i n exper i ments b y 
Brown et al\ (|2005l) : see also iNiemelafc Srecn ivasan (^20_06a) and Ah lers et all (|2009q[ ) . 

The role of boundary thermal properties is also receiving increasing attention in the 
geophysical community in the context of heat transport due to mantle convection. The 
ocean floor and continents impose different thermal conditions at the upper boundary of 
the Earth's mantle: the oceans are well-described as enforcing a flxed temperature, while 
continents act as (partial) insulators, and are modelled as lids of flnite conductivity fully 
or partially covering the convecting fluid. The pre sence of continents is understood to 
affect the convective flow (jGuillou fc Jaupart (1995)), and the effect of finitely conducting 
continents on heat transpor t in mantle c onvection has been investig ated through models 
and numerical simulations ( Lenardic fc Moresi (2003i): [Crigne et~al. (2007a, 6)). 

Careful numerical investigations permit control of extraneous variables that may play a 
role experimentally, thus making it possible to isolate the effect of the thermal boundary 
conditions. Two groups have recently explored this independently : the two-dimensional, 
horizontally periodic computations of I Johnston fc Doeriird (2009) studied fixed t e mper- 
atu re and fixed flux BC s both above and below, while FVerzicco fc Sreenivasan ( 2008h 
and Stevens et al. ( 2010f) compared the effects of fixed flux and fixed temperature lower 
horizontal plates in their cylindrical simulations. No differences between the extremes of 
perfectly conducting and insulating boundaries were observed in either case. However, di- 
rect numerical simulations are as yet unable to attain the high Rayleigh numbers achieved 
experimentally or relevant to, for instance, geophysical or astrophysical applications. 

Analytical upper bounds on convective heat transport: 

In the investigation of transport and scaling properties, mathematical results system- 
atically derived from the differential equations governing the system can play a role. 
The details of turbulent dynamics are beyond the reach of analysis, but bounds on 
averaged quantities can often be obtained, and provide constraints against which phe- 
nomenological theories can be tested, and which are in many situations (though not so 
far in flnite Prandtl number convection) remarkably close to experimental observation. 
In the case of Rayleigh-Benard convection with flxed temperature BCs, a bound of the 

1 /'? 

form Nu < CoRa ' has bee n shown , initia lly with the aid of some plausible statistical 
assumptions (,Ho-^ardl (|l963h : IBuss^ ^M))- More recently, the "background method" 
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introduced in the cont ext of shear flow bv lDoering fc ConstantinI ( 1992 ). motivated by a 
decomposition due to Hopj ( 1941 ). has enable d the above p = 1/2 bound to be proved 
rigorously without any additional assumptions (jDoering fc ConstantinI (|l996l )). This ap- 
proach has turned out to be remarkably fruitful; its applicat i ons t o convection have 



included, among others , studies of porous medium (lOtero et ali (120041) ) . infinit e Prandtl 



number ( Doering et and double diffusive (iBalmforth et "dj (120061) ') convec- 

tion. 

The first investigation to consider the effects of thermal boun dary conditions on rigor- 



ous variational bounds on convective heat transport was that of lOtero et ali (|2002l ). who 



considered upper and lower fixed fiux BCs. This work established an overall bound of 

1 /2 ■ 

the form Nu < CqcRci ' , with the same scaling as in the fixed temperature case; but the 
mathematical structure of the bounding calculations and the intermediate scaling results 
in the two cases turned out to be quite different. When the temperature drop across the 
fluid is fixed, the Rayleigh number Ra is the control parameter, and one obtains bounds 
on the heat transport by controlling the ayeraged heat flux through the fluid boundaries 
from above ( Doering fc Constantin (1996); Kerswell ( 2001 1)). On the other hand, given 
a fixed boundary heat flux, the control parameter R is defined in terms of this imposed 
flux; in this case the averaged temperature difference between the fluid boundaries (and 
hence the Rayleigh numbe r Ra) must b e estimated (from below) in terms of R to find 
bounds on Nu. One finds (|Otero et al\ (2002)) that Nu < ciR^/^, Ra > cai?^/^, unlike 
in the fixed temperature case for which Nu < CqR^^^, Ra = R. It is thus natural to 
wonder how these two extreme cases, corresponding respectively to the idealizations of 
perfectly conducting and insulating plates, are related vis-a-vis their bounding problems, 
and which is more relevant to real, finitely conducting boundaries. 



Outline of this paper: 

In the present work we reconsider the effect of general thermal BCs on systematically 
derived ana l ytical bounds on thermal convection, continuing the program initiated by 



Otero et ali ()2002[ ): we assume for simplicity only identical thermal properties at the top 



and bottom fluid boundaries in the mathematically idealized horizontally periodic case. 

We model imperfectly conducting plates in two different ways. One method is to assume 
mixed (Robin) thermal BCs of "Newton's Law of Heating" type, with a fixed Biot number 
77, and to develop the analysis in a manner which interpolates smoothly between the 
fixed temperature (Dirichlet: rj — 0) and fixed flux (Neumann: rj — oo) extremes (to 
our knowledge the only prior bounding study with general Biot number is the horizontal 
convection work of ISiggers et al\ ( 2004[) . with mixed BCs at the lower boundary). The 
other approach is to consider the more realistic case of a fluid in thermal contact above 
and below with finite conducting plates, restricting ourselves to homogeneous, isotropic 
plates with fixed temperatures imposed at the top and bottom of the entire system. 

In § [3] we formulate the governing equations for Rayleigh-Benard convection and dis- 
cuss various thermal BCs, paying particular attention to the choice of nondimensional- 
ization. Global identities and averages, including energy identities for convection with 
plates, are discussed in § [3l while a bounding principle using the Constantin-Doering- 
Hopf "background field" variational method is derived in § H) The use of a piecewise 
linear background temperature profile and of conservative estimates in § [5] permits the 
derivation of explicit analytical bounds on the Nu-Ra relationship, asymptotically valid 
as i?a cx), as discussed in § |6l For clarity, l3.2H5?T] of the main text treat the case of 
convection with plates, while the corresponding calculations for fixed Biot number BCs 
are presented in a parallel fashion in Appendix |BJ 
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Summary of results: 

For convection with plates, we find that the heat transport depends on d, the ratio of 
plate to fluid thickness, and A, the conductivity ratio, only via the combination a = d/X; 
and that the (conservative) bounding problems with plates and with fixed Biot number 
T] map onto each other when a = rj; this gives a systematic correspondence between the 
"full" problem of conducting plates and the flxed Biot number approximation, without 
stationarity, fixed horizontal wave number or other modelling assumptions. 

Since in general the boundary temperatures are unknown a priori, one must identify a 
temperature scale O extracted from the thermal BCs; a control parameter R, defined like 
a Raylcigh number but in terms of 0, may then be introduced as a measure of the applied 
driving. For sufficiently small Biot number rj (or, equivalently, a), we show that for small 
R we have Nu < 0(R^^'^), Ra > 0{R) as in the fixed temperature case, but that for 
R (and hence Ra) beyond some critical parameter which we estimate as i?t = 0{r)~^), 
we find Nu < ci{r])R^^^, Ra > C2{ri)R^/^ , implying Nu < CnRa^^"^ with intermediate 
scaling as in the fixed flux case. Interestingly, for each 77 > we find = Coo- at least 
at the level of our estimates, the asymptotic scaling in each case is as for fixed flux BCs, 
while fixed temperature BCs give a singular limit of the general asymptotic bounding 
problem. 

Interpreted in terms of convection with plates, the analytical bounds on the Nu- 
Ra relationship confirm that for relatively small R most of the temperature drop oc- 
curs across the fluid. However, for each a > 0, asymptotically as i? — )■ 00 we have 
Nu < c{(t)R^'^ < C^Ra^''^: the bounds scale as in the fixed flux case, providing rigor- 
ous support for the intuition that for large Ra, plates of arbitrary flnite thickness and 
conductivity act essentially as insulators. The asymptotic result Nu < c{a)R^/^, where 
c(cr) = 0{a~^/^{l + 2cr)^/'^) is of particiflar interest, since in this case R may be inter- 
preted as a Rayleigh number in terms of the full temperature difference across the entire 
system,. 



2. Governing equations and thermal boundary conditions 

2.1. Governing differential equations and nondimensionalization 

We consider a fluid of depth h, kinematic viscosity Vf and thermal diffusivity Kf, with 

density pj at some reference temperature To; we also let a be the thermal expansion 
coefficient, Cpj be the specific heat and hence A/ = pf Cpj Kf be the thermal conductivity 
of the fiuid. 

The (dimensional) partial differential equations (PDEs) of motion in the Boussinesq 
approximation, describing the evolution of the fiuid velocity field u* and temperature 
field T*, are 

8n* 1 

— -h u* • V*u* + -V*P* = ufV*^u* + ag{T* - Tq) , (2.1) 
Cr pf 

V* • u* = , (2.2) 

BT* 

— + u* • V*T* = KfV*'^T* (2.3) 

(where g is the gravitational acceleration). In this formulation, the compressibility of 
the fluid is neglected everywhere except in the buoyancy force term, and the pressure 
P* is determined via the divergence- free condition on u*. Variables with an asterisk 
are dimensional, and we take periodic boundary conditions in the horizontal directions. 
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with periods L* and L*, respectively. In the vertical direction, the fluid satisfies no-shp 
velocity boundary conditions u* = at z* = and z* = h. 

The nondimensionalization is chosen to treat the different thermal boundary conditions 
(BCs) at the interfaces between the fluid and the plates at z* = 0, h consistently and in a 
single formulation. For now, we thus let be a general temperature scale, and introduce 
a reference ("zero") temperature Tiot; for given thermal BCs, the approach which turns 
out to be successful is to define O and Tj-cf so that the stationary, horizontally uniform 
perfectly conducting state in the fluid (u* =0, V*T* = C for some constant C < 0) 
takes the nondimcnsional form 

u = 0, T=l-z (0<z<l). (2.4) 

We nondimensionalize using Trof and the temperature scale Q, and with respect to the 
fluid layer thickness h and thermal diffusivity time h'^ / K,f \ that is, we take h, /i^/kj, U = 
Kf /h and pfU^ as our appropriate length, time, velocity and pressure scales respectively. 
For Trcf 7^ To, the nondimensional fluid momentum equation will contain a constant term 
proportional to Trcf — Tq in the direction, which we absorb into the rescaled pressure. 
In summary, the nondimensional variables (without asterisks) are defined by: 

X* . t* u* T* - r,ef 1 p „ r,ef-ro 

x=^, t=- , u=— , T= , p^- 77^ p; z, 2.5) 

h iscal U 6 Pr pfU^ 9 

where igcai = ^I'^f-, U = Kf/h, and Pr and R are defined below. The dimensionless 
periodicity lengths in the transverse directions are = L*/h and Ly = L*/h, and 
A — LxLy is the nondimensional area of the plates. 

The equations for the nondimensional fluid velocity u = {u, v, w) and fluid temperature 
T are thus 

Pr^i (|^ + u- Vu) + Vp = V^u + i^Te^, (2.6) 



dT 
'dt 



V • u = 0, (2.7) 

u-vr = v2r, (2.8) 



with no-slip BCs ujz^o.i = Oj a-nd L^, Ly-periodic BCs in the horizontal x and y direc- 
tions in all variables. 

Here the dimensionless constants are the usual Prandtl number Pr = Vf/nf and the 
control parameter R, defined in terms of the (as yet unspecified) temperature scale Q as 

R^^e. (2.9) 

2.2. Thermal boundary conditions imposed at interfaces 

The specification of the governing equations is completed once conditions on the temper- 
ature at the fiuid-plate interfaces z* — and z* = h are specified. We shall consider both 
thermal BCs applied directly at these interfaces, as in figure 1, and (in ? 12.31 below) the 
case of solid plates in thermal contact with the fiuid; in each case we restrict ourselves 
to fluids with thermally identical upper and lower boundaries. 

Fixed temperature (Dirichlet) conditions: 

The usual and most-studied assumption regarding thermal boundary conditions at the 
interfaces is that the temperature is fixed at the upper and lower fluid boundaries: 

T*U,^o = T*ix*,y\0,t*) = T^, T%,=h = T;. (2.10) 



Bounds for convection with conducting plates 

Upper thermal BCs 



z' = h 




Fluid 



Lower thermal BCs 

Figure 1 . Geometry of Rayleigh-Benard convection system with thermal boundary conditions 
imposed at upper and lower limits of fluid layer. 



These Dirichlet BCs imply a natural choice of reference temperature T^^i = ^T, while the 
imposed temperature drop AT* = — T*|^.^g = T^—T^ introduces a natural temperature 
scale Q — AT*. The nondimensional thermal BCs thus take the well-known form 



T = 1 on X 0, T = on ;z = 1. 



(2.11) 



Fixed flux (Neumann) conditions: 

At the opposite extreme is the fixed flux assumption that the thermal heat flux 
—XfT*, = —XfdT*/dz* through the fluid boimdaries is a constant $. This corresponds 
to the Neumann BCs of fixed normal temperature gradient —f3* at the interfaces: 



T* L.=n = T* 



/ 



(2.12) 



the corresponding temperature scale is = h/S* = h^/Xf, while in this case T^ei is 
arbitrary. In this limit, the dimensionlcss thermal BCs are 



Tr = — 1 on z = and z — 1. 



(2.13) 



Fixed Biot number (Robin) conditions: 

General linear thermal conditions at the boundary of a fluid as in flgure[l]are of mixed 
(Robin) type; in dimensional terms, we write the mixed BCs in the form 

T* + 77*n • V*T* = A;* on z* = 0, T* + 7j*n ■ V*T* = on z* ^ h. (2.14) 

for some given constant < r]* < oo^ These conditions may be interpreted as Newton's 
Law of Cooling (Heating), in which the boundary heat flux is assumed proportional to 
the temperature change across the boundary: -A/n • V*T* = A/(T* - A;*)/?7*. 

We use n — — e^, on z* = 0,h respectively, and nondimensionalize by substituting 
z* = hz, T* = Ticf + 6T. Defining the Biot number rj = r]* /h, we fincjjl 



t The limit rj* 
n-'V*T* +T*/ri* 



e 



0, 



e 



1. 



(2.15) 



oo is treated by writing (|2.14l) in the equivalent form (for i]* > 0) 



0, h, where (for < r;* < oo) = Ai^^/r/* . 



J There appears to be little consensus in the literature as to whether the term "Biot number" 
refers to rj as defined in (|2.15p . or to its inverse ri~^. 
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Figure 2. Geometry of Rayleigh-Benard convection system with conductive plates. 



The so far unspecified reference temperature Ti-ct and temperature scale Q are now deter- 
mined by the condition (j2.4[) on the nondimensional form of the conduction temperature 
profile: requiring T = 1 — z to satisfy the BCs (|2.15p . we find that (for rj < oo) 

e^^l^, T„,^^l^±^^(^t±^. (2.16) 

Having finally fixed a choice of dimensionless variables, the nondimensional mixed ther- 
mal boundary conditions (fixed Biot number) are 

T -riT^ = I + r] on z = 0, T + = -t; on z = 1. (2.17) 

Note that the mixed (Robin) BCs (|2.17p reduce to the fixed temperature (Dirichlet) 
BCs (|2TT|) in the limit ?? -> 0, and to the fixed flux (Neumann) BCs ((2?T3l) in the limit 

— >■ oo; thus we denote 77 = and 77 = cxd as the "fixed temperature" and "fixed fiux" 
cases, respectively. 

2.3. Fluid bounded by conducting plates 

The specification of thermal conditions directly at the fiuid boundaries z* = 0, h, as 
in § 12.21 is an approximation to the experimentally more realistic situation of a fluid 
bounded above and below by conducting plates, with thermal BCs imposed on the plates. 
We consider only the simplest case of plates with equal thickness and thermal properties. 

Beginning with a fluid with properties as in ^ \2.1\ we thus place identical homogeneous, 
isotropic solid plates of thickness hs, thermal diffusivity and thermal conductivity 
As = ps Cp^s Kg above and below the fluid; see figure [5] The spatial coordinates are chosen 
so that 2* = is at the lower boundary of the fluid, and thus the lower and upper plates 
extend from z* — —hs to z* — 0, and from z* = h to z* = h + hs, respectively. 

The governing PDEs in the fluid in the Boussinesq approximation, valid in the region 
< z* < h, are as in (|2.1[) - (|2.3p above, where T* = Tj is the fluid temperature field, 
and the fluid velocity u* satisfies the usual no-slip boundary conditions at z* = and 
z* = h, the interfaces between the fluid and the plates. 
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These equations are coupled to the heat equation for the temperature T* in the plates, 

dT* 

= k,V*^T;, (2.18) 

valid in the lower plate for —h^ < z* < and in the upper plate for h < z* < h + hg. 

At the interfaces (where n — e^) we require continuity of temperature T* and normal 
heat flux An • VT* = XdT*/dz*. With T*, T*^ and T*,^ representing the (dimensional) 
temperature in the fluid, lower plate and upper plate, respectively (letting subscripts I 
and u identify the plates), continuity of temperature may be written as 

p,l\z*=Q — J-f\z'=0 , J-f\z*=h — -'-p,u\z*=h , 

and similarly for flux continuity. However, it is more convenient to treat T* as a single 
temperature field, continuous but with discontinuous derivative, which coincides with T* ^ 
for —hs < z* < 0, with for < z* < h, and with T^* „ for h < z* < h+hs] and we write, 
for instance, T*|,.=o+ = lim,.^o+T* = T}l,=o.oi {dT* / dz*)l.=^_ = {dT;jdz*%>=o 
(see Appendix VK\ concerning notation). We may then express the continuity of tempera- 
ture and heat flux at the fluid- plate interfaces as: for each x* , y* and t* , 



and 



^ dT* 

* dz* 



T*U.=o-=T*U.=o+ , T*U.=h-=T*U.^h+ , (2.19) 

dT* 



z*=0- 



dT* 



z*=0+ 



dT* 



z'=h- c)z* 



(2.20) 



z*=h+ 



We assume that the entire plate-fluid-plate system has Dirichlet thermal boundary 
conditions in the vertical direction (in addition to horizontal periodicity in all variables) , 
with fixed temperatures at the bottom of the lower plate and the top of the upper plate, 

r*U.=-,,, = r;, t*|,.=^+,,^ = t: ; (2.21) 

and we define the overall temperature drop across the system as 

A*=T;*-T*. (2.22) 

Nondimensionalization: 

The coupled governing PDEs are nondimensionalized with respect to the fluid param- 
eters, as described previously in (|2.5p . As before, the rescaling yields the dimensionless 
Prandtl number Pr = i^f/nf, and the parameter R defined as in (|2.9I) . This R will be 
our control parameter, in lieu of the usual Rayleigh number, because the latter is defined 
in terms of the temperature drop across the fluid, whereas a priori we know only the 
temperature drop A* across the entire system (|2.22[) . 

The presence of the plates introduces as additional parameters the nondimensional 
plate thickness, thermal diffusivity and thermal conductivity — equivalently, the plate- 
to-fluid thickness, diffusivity and conductivity ratios — 

d=^, n=^, A=^; (2.23) 

h Kf Xf 

we also have the density and specific heat ratios p — p^/pf, Cp — Cp^s/cpj, where pcp = 
X/k. We now introduce the ratio a of the dimensionless thickness and conductivity, 

d hs Xf „ .N 

this will turn out to be the main physical parameter of the problem, playing an analogous 
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role to the Biot number 77 of ()2.17p Lastly, we need to choose the reference temperature 
Trcf and temperature scale 0, so that in nondimensional form the temperature field is 
T — (T* — Ti.of)/6; the imposed boundary temperatures (j2.2ip then become 

Tu = TU=i+d - ^^^Q^^ Ti = TU_, = -^-^ = ^ + (2.25) 

It is again convenient and consistent to define Q and Tjcf so that the dimensionless linear 
conducting state in the fluid (0 < z < 1) is given by (|2.4I) . By flux continuity (see (|2.3ip ). 
we have — — 1/A in the plates, so that the dimensionless temperatures at the lower 
and upper boundaries of the system are T/ = T\z=a — d{—^/^) = 1 + cr, ?u = —a, with 
total overall temperature drop Ti-Tu = {T* -T*)/(d = l + 2a. Substituting into (P?^ 
and solving for T^cf and 0, we conclude that appropriate choices are 

A* T* ~ T: T* + a{T*+T:) 

° = TT^ = TTmTa' ^ ■ ^ ^ ^ 

Dimensionless formulation of Boussinesq convection with plates: 

The nondimensional formulation of the governing PDEs and BCs for Rayleigh-Benard 
convection with conducting plates is now complete: The equations for the dimensionless 
fluid velocity u = {u,v,w) and temperature T = T/, valid on < z < 1, are (|2.6p - 
p.7p with no-slip vertical velocity BCs, exactly as before. The continuous (piecewise 
smooth) temperature field T satisfies an advection-diffusion equation in the fluid, and 
heat equations in the plates, so that we have 

dT 
'dt 

dT „ 

dT 
'dt 

The dimensionless interface and boundary conditions are: at the fluid-plate interfaces, 
we have continuity of temperature 

T\,^o- = TU^o+, = r|,=i+ (2.30) 

and of heat flux 



kV^T, 


-d < z <0 


(r = 




(2.27) 


V^T, 


< z < 1 


(T = 


Tf), 


(2.28) 


kV^T, 


1 < z < 1 -1- d 


(r = 


T ) 


(2.29) 



A — 

dz 



dT 



dT 
'dl 



- A — 
dz 



(2.31) 

1+ 



=0+ 

while the applied temperatures at the upper and lower boundaries of the system are 

TU=^d = Ti = 1 + a, T|,=i+<j T„ = -a. (2.32) 

In proceeding further, the formulation of global identities and of a bounding principle 
for these coupled equations in the plates and fluid is greatly simplified by an appropriate 
well-chosen notation; we relegate some of our notational definitions to Appendix [X] 

f Not e that a = d/ X is sometimes referred to as " the Biot number" of a system; see for 
instance ISparrow et al\ (llQG^ ): IChapman et al\ l|l98Cf ): iGrigne et al\ (|2007af ). However, we use 
the term Biot number specifically to denote the constant ri in given (mixed) thermal BCs of 
the form T + rjn ■ VT = const, applied at the fluid boundaries. When plates are present the 
Biot number then depends on a perturbation horizontal wave number; and d/X is in fact the 
Biot number at zero wave number, or that appropriate to the thin-plate limit; see for instance 
ifCross fc HoficnbcrE.199A Section VIII. F.l) (and also the previous footnote). 
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Limiting values of a: 

It is instructive to consider the interpretation of the hmits cr — > and cr — oo, when 
(for fixed fluid height h and conductivity Xj) either the plate thickness hg or conductivity 
Xs approach or cx)l|| 

In the limit of vanishing plate thickness /i^ — > 0, according to (|2.2ip the temperatures 
are fixed at the lower and upper boundaries of the fluid. Similarly, when the plates are 
perfect conductors, As — ?> oo, they sustain no temperature gradient, and the temperatures 
at the fluid boundaries coincide with those applied to the plates. In both of these cases, 
d — >■ and A — >■ oo, we recover the fixed temperature BCs ()2.1ip . so that a — d/X — >■ 
corresponds to the fixed temperature limit; the corresponding temperature scale is just 
given by the applied temperature drop, 6 = limcr^o A*/(l + 2a) = A*, as expected. 

Somewhat more care is required for a — > oo, as by (I2.26P we then simultaneously 
need A* = Tj* — T* oo for 9 to remain finite. Since then 8 — limo— i.oo A*/(l + 2a) = 
limo-^^oo (A*/2(t), this implies that lim^-^oo {XgA* /2hs) — XfQ/h is finite, while A* — >• oo 
and either As ^ or /i^ — >■ oo. Thus $ = limg-^oo (AsA*/2ft,s) is well-defined, and is the 
magnitude of the fixed imposed flux across the system. 

Specifically, the limit As — )■ corresponds to perfectly insulating plates; in this case, 
by (|2.20l) the vertical temperature gradient T*. across the plates must diverge so that 
XsT*,\z'=o- = XfT*,\z'=o+ remains bounded, and equals the boundary flux — $ (sim- 
ilarly at z* = h). A lternatively, for < As < oo, we may model infinitely thick plates 
( Hurle et Ml (119671) ) l)y letting hg oo and A* — >■ oo so that the global temperature 
gradient limft,^^oo (^A*/(/i + 2hs)) = lini/j^^oc A*/2ft,s) remains finite, and hence so 
does the overall flux lim/i^_>oo (— AsA*/2/is) = — $. In either case A — >■ or d — >■ oo, we 
have a — d/X — >■ oo, which gives the fixed flux limit with BCs (j2.13p : and the temperature 
scale is chosen as 9 = h^/Xf. 

The limiting cases — > and u — > oo are thus best treated by im posing the thermal 
BCs on the fluid boundaries as in ii l2.21 as in the literature (for instance lDoering fc Constantin 
(l99fih : lOtero et al\ (j2002f )). In the following we consider plates of finite thickness and 
conductivity, so that < cr < oo, and (|2?6 l) - ([2Jl) and (|2?27l) - (|232|) apply. 



3. Global identities 

We next derive some exact relations between averaged quantities, using the notation 
outlined in Appendix|3 First we need to recall the definitions of the Rayleigh and Nusselt 
numbers, as the relationship between these is the primary goal of our investigation. 

3.1. Rayleigh and Nusselt numbers 

Rayleigh number: 

We define the nondimensional horizontally- and time-averaged temperature drop across 
the fluid as 

= -{T\l^o) = (rU=o - TU.i) ^ (3.1) 

where AT* = (T*|^*=o ^ T*]^-^^}/) (this is well-defined in the presence of plates since T 
is continuous at the interfaces (12.301) ). We observe that this temperature difference AT 
is known a priori only for fixed temperature BCs (or equivalently, when 77 = or cr = 0), 
in which case AT* = 9, AT = 1. The conventional Rayleigh number Ra is defined in 

t We do not consider situations where ha and As approach and/or 00 simultaneously. 
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terms of the averaged fluid temperature drop AT* as 

Ra ^AT* = ^AT, (3.2) 

and is related to the control parameter R (defined in (j2.9p in terms of <d) by 

Ra^RAT. (3.3) 

Nusselt number: 

The Nusselt number Nu is a nondimensional measure of the enhanced vertical heat 
transport across the fluid due to convection, relative to the conductive heat transport 
associated with the same temperature drop AT* . Its expression in terms of flow quantities 
is standard: one writes the thermal advection equation in the fluid (|2.8p as a conservation 
law, Tt + V • J = (using (|2.7p ). where the dimensionless heat current J = Jc + J« is 
the sum of the conductive and convective heat currents, Jc — — VT and Jy = uT. Then 
Nu is defined as the ratio of the total (averaged) vertical heat transport, (/^ Sz - J), to 
the purely conductive transport (/^ • Jc), to give the well-known expression 

A more useful formula, which allows us to estimate Nu from the equations of motion, 
is found by relating {Jj wT) to the time-averaged temperature drop and boundary fiux. 
To do so, we begin by taking the horizontal average of the temperature equation ()2.8|) . 
using the horizontally periodic BCs, to get 



Tt + V ■3 = Tt + {wT -T,)^=Q. (3.5) 
Integrating over z and using the vertical no-slip boundary conditions on w, 

^jT + A[\^-T.)^dz^^jT + A (-T.) = 0. (3.6) 



dtjf Jq ^^Jf 



Now one may show, using techniques similar to those i ntroduced b y [l5oering fc Constantin 



(|l992| ) in the context of shear fiow (based on an idea of iHopj (119411) ). i;hat the fiuid ther- 
mal energy ||T|| ^ = is uniformly bounded in time; for Rayleigh-Benard convection 



with fixed temperature BCs this boundedness was verified bv lKerswelll (|200ll) . It follows 
via JjT < A^l"^ 11-^11/ ^^^^ // ^^'^ uniformly bounded. Hence on taking a time aver- 
age of p.6p . the time derivative term vanishes, and we find (,—T^^_^_^^ = 0, expressing 
the expected result that, on average, there is a balance between the heat fluxes entering 
the fluid layer at the bottom and leaving it at the top. 

This motivates the definition of /3, the nondimensional horizontally- and time-averaged 
vertical temperature gradient, or equivalently, the nondimensional heat flux, at the in- 
terface between the fluid and the plates: we define 

Note that this quantity is known a priori only for fixed fiux BCs (or equivalently, in the 
limits 77 = oo or cr = od), in which case (3 — 1. In the presence of plates, by (j2.3ip we 
also have 

/3 = A(-T,)|^^„_-A(-T,)|^^^^. (3.8) 
If we had a general uniform bound on T, we could immediately take a time average 
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of p.Sp and deduce that (Tt) = 0. However, for fixed flux BCs we have n o maximum 
principle on T to provide such an a priori bound. Instead, foUowing Otero et al. (.2 0021 ). 
uniformly in thermal BCs we multiply (13.51) by z and integrate to obtain 



I zTdz+ j z{wT -T^)^ dz = ; (3.9) 

and as before, via /g zTdz < (3A)-i/2 ||r||^ and the uniform boundedness of ||T||p the 
time average of the first term in p.9p vanishes. By integration by parts and the no-slip 
BCs, the second term in (|3.9p becomes z [wT — T^) ^ dz — —Tz\z=i- — -j Jf wT + 
T|^_i^; taking time averages of (|3.9p and using p.ip and p.7p . we obtain 

1 
A 



wT } ^ P - AT. (3.10) 



/ 

Substituting p.lOp into p.4p , we now obtain the fundamental Nusselt number identity, 

Nu = 4-; (3.11) 
AT ^ ' 

while via p.3p . Nu, Ra, and the control parameter R are related by 

NuRa^Rp. (3.12) 



Organizational remark — Biot number calculations in AvvendixW[ 

In the following sections we extend the bounding principle, previously studied in the 
fixed temperature and fixed flux extremes, to our more general thermal boundary condi- 
tions. As described in §§[12H131 we model imperfectly conducting fluid boundaries in two 
ways: by imposing mixed BCs of finite Biot number, and by assuming the fluid to be in 
thermal contact with (identical) plates of finite thickness and conductivity. Since details 
of the calculations differ in these two cases, for clarity of presentation we have separated 
them: in the following sections of the main text we consider convection with bounding 
plates, while the analogous results for finite Biot number are relegated to Appendix [B] 



3.2. Relation between /3 and AT for convection with plates 

In the general case, when the boundaries of the fluid are neither perfectly conducting 
(fixed temperature) nor perfectly insulating (fixed flux), neither AT nor /3 is known a 
priori. However, they are related via the thermal BCs; this is crucial to formulating a 
bounding principle on the Nusselt number, as once one of /3 and AT is estimated, the 
other and, using (|3.1ip . hence Nu may also be controlled. 

For convection with bounding plates, taking horizontal and time averages of the heat 
equations (|2.27p and p.29p . we find that in each of the two conducting plates 

k{Tzz) = (Tt) = (3.13) 

(using a maximum principle on T for a < oo); consequently the averaged temperature 
gradient (Tz) is a ^-independent constant in each plate, separately for —d < z < 
and 1 < z < 1 -I- d. In particular, in the lower plate this gives {T\z=o — T\z=-d)/d = 
{Tz)\ = (where in the last identity we used (13. 81) ). or 



(TU^o) = {Tl=-d) ~d = Ti~ al3. 



(3.14) 
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Similarly, in the upper plate we find {T\z=i+d — T\z=i)/d ~ {'^z)\^^-^^^, or 

(T|,=i) = (TU^i+d) +/3j=Tu + a/3. (3.15) 

Subtracting p.lSp from p.l4p . and using p.ip and (|2.32p . we obtain the basic relation 
between AT and (3 for conducting plates, 

AT + 2(7/3 = 1 + 2o- (3.16) 

(compare the analogous result (jB 1[) for fixed Biot number). 

3.3. Energy identities 

We next obtain the basic "energy" identities from the governing Boussinesq PDEs, 
which allow us to relate Nu to the momentum and heat dissipation. In evaluating time 
averages, we again use the fact that u and T are a priori bounded in L^. 

Kinetic energy: 

The kinetic energy balance is obtained by taking the inner product of the momentum 
equation (j2.6l) with u; standard integration by parts, using no-slip BCs and incompress- 
ibility, and time averaging yields the identity across the fluid (also using p. 101) ) 



1 

R 



f 



Vup ) = ( / wT)^A{P- AT). (3.17) 



Observe that (|3.17l) implies that (3 > AT, so that by (13. lip we have Nu > 1, as expected. 

In the presence of finitely conducting plates, by p.l6p we can solve for one of AT and 
P and state the energy identities in terms of the other. We shall state our results (for 
fT < oo) in a way that permits the derivation of an upper bound on /3; this formulation, 
suitable for small a, reduces to the known fixed temperature identities as cr —>■ 0. Thus, 
using ((3l6)) in the form /3 - AT = (1 + 2cr)(;3 - 1) to substitute for AT, (|3T7)) gives 

i({|Vup})=yl(l + 2a)(/3-l), (3.18) 

where we have also used the weighted integral (jA 7p . defining u = in the plates. 

Thermal energy: 

The presence of plates modifies the global thermal energy balance, since the thermal 
BCs (|2.32l) are given at the ends of the plates, not of the fiuid. Multiplying (I2.28P by 
T, integrating over the fluid, integrating by parts and taking time averages, we find the 
general thermal energy identity over the fluid, 

^VT\^^=a(Tt:\IZ^^), (3.19) 

using the notation introduced in (|A 31) . Beginning with (I2.27P and (j2.29p and proceeding 
similarly over the plates, we find 

I VTI^ )=A(^ TT;C_,) , I I VT| A =aU T%Ct+) ■ 



We now multiply the identities over the plates by pcp = X/k before adding them to 
the fiuid identity (|3TT9l) : since from (|2:30| and ([OT]) we have ATT^|^=o- = TT^\z=o+ 
and TTz|z=i_ — ATTz|z=i_|_, all terms evaluated at the fluid-plate interfaces cancel by 
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the temperature and flux continuity conditions. Thus we find 

-A{xm\ltt,)- (3.20) 

To evaluate the boundary terms in p.20p . we use the known values of T at 2; = —d 
and 1 + d (I2.32p . and the result from (|3.13l) that the averaged temperature gradient 
(Tz) is constant in each plate; using p.8p we find {Tz)\z=-d = {Tz)\z=o- = and 
(5^^)12=1+^ = -/3/A. Writing the left-hand side of ((3^ using the shorthand (lA 7|) for 
the weighted integral over the entire plate-fluid-plate system, we substitute the boundary 
conditions to obtain the global thermal energy identity 

{{\VT\''})=A/3{Ti-T^) = A{l + 2a)p. (3.21) 



4. Background fields and formulation of bounding principle 

4.1. Background flow decomposition 
The C onstantin-Doering-Hopf "background" method for the convection problem (jPoering fc Constantin 



( 1996[ )) relies upon a decomposition of the temperature field r(x, t) across the entire sys 



tem into a background profile f = f(z) which obeys the inhomogeneous thermal bound- 
ary conditions, and a space- and time-dependent component 0(x, t) with homogeneous 
boundary conditions: 

T{x,t) ^f{z) + 9{x,t). (4.1) 
For the v elocity decompo sition, the assumption of zero background flow is likely to be 



sufficient (jKerswelll ([20011)). It can nevertheless be helpful to introduce a "fiuctuating" 
field V over which we shall optimize, conceptually distinct from the velocity field u solving 
the Boussinesq equations; so we write u(x, t) = v(x, t) = {u, v, w). 

The function f (z) is for now arbitrary, provided it satisfies the boundary and interface 
conditions on T; that is, from (|2.30p - (|2.32|) we require 

f{-d) ^Ti^l + a, f(l + d) = T„ = -a, (4.2) 

and 

f(0-)=f(0+), Af'(O-) = f'(0+), f(l-)=f(l+), f'(l-) = Af'(l+) (4.3) 

(note that if A 7^ 1, f (z) has discontinuous slope at the fiuid-plate interfaces). When the 
upper and lower plates are identical, it is sufficient to consider only symmetric background 
fields satisfying f'(0+) = f'(l— ) (compare p.Tp ): we define 

Af = f(0)-f(l), 7 = -f'(04-) = -f'(l-), (4.4) 

and observe that by (|4.3p we have — f'(0— ) = — f'(H-) = 7/A. 

Since the background f carries the same boundary and interface conditions as the 
temperature field T, the fiuctuation 9 = T — f vanishes at the outer ends of the plates, 

ei^-d = 0\z=i+d = 0, (4.5) 
and also satisfies the temperature and flux continuity interface conditions, 

0\z=0- — &\z=0+, ^&z\z=0- — (^z\z=0+, Qz = l- = Q\z=l+, dz\z=l- — ^&z\z=l+- (4.6) 



Substituting the decomposition T 



— T + 9 into the Boussinesq equations with plates 
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(IM])-(I2I7]), (|2:27| - ([2:29l) . we obtain the PDEs for the fluctuating fields: 



Pr ^ ( ^ + ■ 



86 



Vp = \7^v + R9e^, 



V • V = 0, 

89 



KT 



vV0 

89 



< z < 1 , 


(4.7) 


< z < 1 , 


(4.8) 


- d < z < , 


(4.9) 


< z < 1 , 


(4.10) 


1 < 2 < 1 + d . 


(4.11) 



where in (j4.7l) we have absorbed the Rf term into a redefinition of the pressure p. 
Here v satisfies no-shp BCs v|^=o,i = and can be defined across the entire domain 
z G [—d, 1 + d] by setting v = in the plates. 

4.2. Energy identities for fluctuating fields 

The evolution equation for the field 9 is obtained in a similar way to (|3.20p : multiplying 
each of (|4.9p ~ (|4.1ip by 9, integrating over the relevant domains, integrating by parts 
(using (|4.8|) ). multiplying the integrals over the plates by pcp = X/k and adding the 
results, we find 



1 d 
2dt 



A / \V9\ 



\V9\ 



A 



99, 



\z=0+ 







- 9,f' + xj 9,f' 









A 



At 



I z— — d 



I 2=0- 



+ Af'' 



I l+d 
\z = l+ 



771 l + d 



-{\V9\^]-{9,f'}- jw9f'. 



w9f' 
(4.12) 



No boundary terms remain, since all the terms at the fluid-plate interfaces cancel due to 
the continuity conditions (j4.3|) and (14.61) . while the extremal boundary terms A ^^z|^__^ 

and A t'6'|^^^^ vanish by the homogeneous Dirichlet conditions (|4.5I) on 9. 

An identity between the norms of gradients of T and 9 will permit us to relate the 
fluctuating field 9 to the unknown fiux /3 (via (|3.2ip ): the decomposition (|4.ip implies 
|VTp = IV^+e^f'p = \V9\'^ + 29^f'+f'^, and taking the conductivity-weighted integral 
(|A 7[) we obtain 

{|VT|2} = {\V9\^} +2{9,f'} + {f'2} . (4.13) 

We eliminate the {9zf'} term by adding 2 • (I4.12p + (|4.13p : time averaging, we find 

({IVTp}) = - {{\V9\'}) - 2 {{f'w9}) + {f'2} . (4.14) 

The relation u = v between the velocity field u and fluctuations v is incorporated into 
the upper bounding principle in the form 



^({|Vup}) = i({|VvP}) 



(4.15) 
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Balance parameter and quadratic form: 

In order to formulate upper bounding principles for the Nusselt number, we now take 
ap propriate linear combin ations of the above identities, using a "balance parameter" 
b (jNicodemus et all ( 1997 )). (When the evolution of the norm of v is also taken into 
account, in general such linear combinations may in fact contain up to three free param- 
eters, over whic h one might op t imize to obtain the best possible bound available within 
this formalism (jKerswelll (|l997l . 120011) ): we shall not pursue this generalization here.) 

Forming the linear combination b ■ (|4.14p + {1 — b) ■ (|4.15p . we obtain 



6({|vrp}) 



1-6 
R 



({|Vu|2}) = fe{f'2} - (b{\ve\^} + 2b{f'we} 

^b{f'^} ~bQr,RA^,0], 



b-i 

R 



{IVvp} 



(4.16) 



where we define the quadratic form in the presence of plates, using a weighted integral 
across the system, as 



Qf [v, I 



1 



bR 
V6i|2 



Vvl- 



2f'we 



^|Vv| 



2 + 2f'we 



\ve\ 



n jf 

Here we have defined an "effective control parameter" Rc via 

b 



+ \ |V 



6-1 



(4.17) 



(4.18) 



having observed that Qf,R^ depends on R and 6 only through the combination bR/{b— 1). 
We desire a positive balance parameter 6 so that a lower bound on Qf.R^ should imply 
an upper bound on /? and/or a lower bound on AT; since i?c > is necessary for Q.f,R„ 
to be a positive definite quadratic form, we thus require 6 > 1. 

Now substituting the identities p.lSp and (I3.2ip for the momentum and thermal dis- 
sipation into (I4.16|) . we obtain after rearranging 



(l-t-2a)(/3-l)-6(^-{f'^}-(l + 2a)j --Q,.fljv,^?] (4.19) 
(we could alternatively substitute (3—1 = (1 — AT) /2a to get an expression only in AT). 



4.3. Admissible backgrounds and a bounding principle 

Although the relation (|4.19l) is exact, it does not permit us to compute (3 since we do 
not have access to sufficient analytical information about the fields v(x, t) and 6{x, t) 
solving (|4.7|) - (|4.11|) . The basic idea of the background fiow method for obtaining upper 
bounds is that, given R, if for some f and 6, Qf,R„ can be shown to be bounded below, 
then (|4.19p yields an upper bound on and ultimately (v i a p. 161) and p. lip ) an upper 
bound on the Nusselt number Nu ( Doering fc ConstantinI ( 19961) ). 

Furthermore, by widening the class of fields v, 9 over which the minimization of Qf,R^ 
takes place (provided this class contains all solutions of (|4.7p - (|4.1ip ). a (weakened) lower 
bound on Q.f,R,, (which, if it exists, must be zero) may indeed be demonstrated. Note that 
if the dynamical constraints on v and 9 imposed by the governing PDEs are removed, 
so that no assumptions are made on the temporal structure of these fields, it is sufficient 
to minimize Qf,R^ over stationary fields. We thus consider, and denote as allowed fields, 
scalar fields 9{x) and divergence-free vector fields v(x) which satisfy the (homogeneous) 
boundary and interface conditions consistent with the given problem; in our case of 
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convection with plates, these are horizontal periodicity for v and 6, the no-slip condition 
V = at z = 0, 1, and that 9 satisfies (|431) - (j46l) . 



Admissible and strongly admissible backgrounds: 

For each > (that is, for each R > and b > 1), we call a background field f(z) 
admissible if it satisfies the same boundary and interface conditions as T, in this case 
(|4.2[) - (j4.3p : and if the resultant quadratic form Qf,R^ is nonnegative, Qr./jjv, 6'] > for 
all allowed fields v and 9. 

Consider now again the quadratic form Q.f.RA'^^^]^ which from (I4.17P may be written 
(for stationary fields) as 

&A,[v,0] = Q^,fljv,^?]+A [\V9\^ + X I \V9\\ (4.20) 

J I J u 

where the quadratic form Qf,R^ is defined as an integral over the fiuid layer only, as 

1 



Qf,RA^,o] 



— iVvp + 2f'we' + ive'p 

Rr, 



(4.21) 



Note that Qt.r^ depends on the background f only through its values on the fluid domain 
< z < 1, that is, only on its restriction r = f |[o.i]- Now since the contributions to Q,f,R^ 
from the plates are clearly nonnegative, from (j4.20p we immediately deduce 

Qf,RA^,e]>QrMA^,9]- (4.22) 

that is, a lower bound on Qf,R^ implies a lower bound on Qf.Ro- 

This motivates the definition of a stronger condition on the background f{z) sufficient 
for obtaining an upper bound, in which we require positivity of the quadratic form over 
the fiuid alone, without assistance from the plate contributions. Correspondingly, we 
enlarge the class of fields over which we minimize: Since we do not have much control 
over 9 at the fluid boundaries, we shall leave the BCs on ^ at z = 0, 1 unspecified. Thus we 
say that f(z) (satisfying the appropriate boundary and interface conditions) is strongly 
admissible if Qt^r^['v,0] > for all sufficiently smooth horizontally periodic fields v(x) 
and 6{x), where v is divergence-free with v = at z = 0, 1. Clearly, by (|4.22p strong 
admissibility implies admissibility. 

Our analysis in § 15.31 below shall in fact yield a condition for strong admissibility on 
the piecewise linear background field fs (or equivalently, on its restriction to [0, 1]), so in 
the following we restrict ourselves to studying this condition. 

Fourier formulation of strong admissibility condition: 

Due to the horizontal periodicity of the problem, we may reformulate the strong admis- 
sibility condition Qf,R^['v, 9] > in horizontally Fourier-transformed variables. To do so, 
we Fourier decompose the vertical component of velocity w = • v and the temperature 
fluctuation 9 in the usual way, 

wix, y, z) = ^ e'('=="^+'=»^)u;k(^), 0{x, y, z) = ^ e'''''-'=+''yy^9k{zy, (4.23) 

k k 

here we use the notation k = {kx,ky) — {'^'^nx/ Lx,2T:ny/ Ly) for the horizontal wave 
vector, with fc^ = |kp; we also write 9^. for the complex conjugate of ^k, and D = 
d/dz. We can use incompressibility to express the transformed horizontal components of 
velocity in terms of the vertical component, so that the admissibility criterion may be 
written completely in terms of the Fourier modes i&k and 0k- This considerably simplifies 
the formulation, particularly since different horizontal Fourier modes decouple in the 
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quadratic form Qf.R^- For strong admissibility we do not impose BCs on 0k at z = 
0,1, while the no-slip bomidary condition and incompressibility imply that the BCs 
on w\^{z) are — DtBk = for z = 0,1. We note also that wq — 0; this follows 
from incompressibility and horizontal periodicity via Awz = Jfj^Wz dxdy 

Vy) dxdy = 0, which implies using w\z=o = that wJ = for all z. 

Substituting ()4.23|) into (|4.21|) and using incompressibility, as in lOtero et al\ we 
can write the quadratic form Qf,R^ evaluated on allowed (stationary) fields v and 9 as 



1 



ivvr 



2f'we + \ve\ 



> 



A^Qk, 



where (see Constantin fc Doerind (1996); Kerswell (2001)) 



1 



•2|Du;k|' 



dz 



(4.24) 

2f'Re[wk^k] 

(4.25) 



note that (|4.24p is an equality for two-dimensional flows. 

Since the class of fields v and considered for strong admissibility includes fields 
containing a single horizontal Fourier mode, it is clear that Qf,R^ is a positive quadratic 
form if and only if all the quadratic forms Qk = Qk;r,i?„ are positive. Thus the strong 
admissibility criterion for background fields f{z) (for given i?e > 0) may be formulated, 
in Fourier space, as the condition that Qk[w'k,6'k] > for all k and for all sufficiently 
smooth (complex- valued) functions Wk(z), 9\^{z) satisfying = Dwk = at z = 0, 1. 

Bounding principle: 

The expression (|4.19|) now implies an upper bounding principle for the Nusselt number: 
for each ii > 0, if we can find a 6 > 1 and an admissible background field t(z) (so that 
Sr,fle['^i^] ^ for all allowed v and 0), then the averaged boundary heat flux /? is 
bounded above according to 



= 1 




' dz 



l+d 



dz + \ f"dz\= Ba [f ; b] , (4.26) 



while the identity (j3.16p then implies a corresponding lower bound on the averaged 
temperature drop across the fluid AT, 

b 2a 



AT > 1 + 2(t6 



{f'2}EEP.[f;&], 



(4.27) 



Al + 2(j 

Via (13. lip , together these bounds yield an upper bound on the Nusselt number: 

Nu<Af„[f-b]^B,[f;b]lV„[f-b]. (4.28) 



5. Piecewise linear background and elementary estimates 

For each R, the best upper bound on the Nusselt number achievable in the formula- 
tion developed above is obtained by optimizing the upper bounds Na[T]b] (14.28^ over 
all admissible backgrounds f{z) and balance parameters & > 1. Careful numerical stud- 
ies obtaining such optimal solutions of analogous bounding problems have been per- 
formed for plane Couette flow (which is relevant to fixed temperature convection) by 
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Figure 3. The piecewise linear background profile ts{z), witli t' = — 7/A in the plates, 
f' = —7 in the fluid boundary layer, and f ' = in the bulk. 



Plasting fc Kerswelll ( 2003 ) and for infinite Prandtl number convection by lerley et al 
(l200d) . 

Ratlier than attempting such a full solution of the optimization problem for the upper 
bound, though, we consider only a restricted class of profiles f(z), for which we shall 
enforce the strong admissibility criterion through Cauchy-Schwarz estimates; we thereby 
much more readily obtain explicit, albeit presumably weakened, analytical upper bounds 
on Nu for Rayleigh-Benard convection with conductive plates. 

5.1. Piecewise linear background profiles in presence of plates 

Following Doering fc ConstantinI ( 19961 ) and subsequent works, we introduce a family of 
continuous, piecewise linear background profiles fs{z) parametrized by (5 (0 < 5 < 1/2), 
for which in the fluid fg — —7 for < z < S and 1 — (5 < z < 1. By the interface conditions 
(|4.3p . the (constant) gradient in the plates is then given by = —7/A for —d < z < 
and 1 < z < 1 + d, so that we define fs as follows: 



t{z) = Ts{z) 



7(5 - 
7(z - 



72/A, 

-s), 



7(5-7(2-l)/A, 



-rf < z < 0, 
0<z<S, 
S < z < 1 - 5, 
1-6 <z <1, 
1 < z < 1 + d, 



(5.1) 



where we still need to find 7 and the average fa in terms of S and the parameters in the 
problem; see figure [3l 

The intuition behind this definition is that in order for f(z) to be strongly admissible, 
the indefinite term 2f' w9 in Sf.fle [v, 0] (see (I4.21|) ') should be controlled by the other, 
positive terms. With this choice of background, 2f'w9 vanishes in the bulk of the fluid 
domain, and is nonzero only near the fluid boundaries, where w and Wz are small. Fur- 
thermore, since f' is piecewise constant, explicit analytical bounds are readily attainable, 
giving (non-optimal) rigorous bounds on the Nusselt number. 

Observe that in the fluid region < z < 1, fs{z) is reminiscent of observed mean 
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temperature profiles in convection, with strong gradients in a narrow thermal boundary 
layer of width ^ (5bl near the boundaries and approximately constant temperature in 
the bulk. This suggests that 5 might be interpreted as modelling the thickness of the 
thermal boundary layer (see also the discussion following (|6.14p below). 

Since the background fg defined in (IS.ip should satisfy the BCs (I4.2p . we must have 
7=5 (-d) = fa + 7(5 4- 7d/A = fa + ^{S + (t) = 1 + and fs{l + d) = Tq - + cr) = -ct; 
solving for 7 and for fa (for cr < cx)), we find 

1 + 2(7 _ 1 , , 

completing the specification of the background fs{z). We can now compute 
Af = fsiO) - fsil) - 27(5 - ^(i±M = 1 + 2(7(1 - 7) 

(5 + (7 

(compare (jB 5[) ): since (5 < 1/2, we remark that 1 < 7 < 1/26 and Af < 1. It follows 
also that fs{0) = fa + 7(5 = 5(1 + Af) = 1 + (7 - (77 and f5(l) = i(l - Af ) = -cr + (77, 
which shows that 

fs — a-fg — 1 + (7 at z = 0+ , f^- + afg — —a at 2; = 1— : (5-3) 

at the fluid boundaries z — 0+ and z = 1 — , the piecewise linear background in the 
presence of plates fs{z) satisfies the mixed thermal BCs (j2.17p with Biot number rj = a. 
To evaluate the bound (|4.26p . (14.271) for this background profile, we compute 

i{f'2}=Ay" f'^dz + J f''Uz + xJ f''Uz = Xd(^^y + 26f + Xd(J^y 

= 2f{6 + a)=j{l + 2a). (5.4) 



Substituting, the upper bound (|4.26l) on (3 and lower bound (|4.27p on AT in the presence 
of conductive plates using a piecewise linear (pwl) background then take the simple form 

1 2(5 

/3 < Bp^iAS, b) = BAfs;b] = l + b^^^^ = 1 + 6(7-1), (5.5) 
AT > Ppwi..((5, b) = t>,[fs;b] = l- b "^] 7 = 1 + 6(Af - 1); (5.6) 

O + (7 

and the corresponding upper bound on the Nusselt number is Nu = /3/ AT < A/'pwi.o-((5, b) — 
'Bpw1,<t((5, 6)/2?pwicr((5, 6). Since 6 > 0, these bounds satisfy Bpvjia{S,b) > 1, I?pw1,ct((5, 6) < 
1, and hence Afpwi,a-{S,b) > 1, as one might expect. Observe that the bounds Bpv,i.a-i6,b) 
and 25p^i £r(i5, 6) do not depend explicitly on the control parameter R, but rather indi- 
rectly through the admissibility condition on 6. It remains, in § 15.31 to find conditions 
on (5 for which fs{z) is (strongly) admissible. 

5.2. Correspondence between bounding problems with and without plates 
The preceding 13 . 2f[?TT] concern the formulation of a bounding principle, and the deriva- 
tion of explicit formulae for the bounds on /3 and AT in the case of piecewise linear 
background fields fs, for Rayleigh-Benard convection in a fluid bounded by conducting 
plates with dimensionless thickness d and conductivity A. In Appendix [Bl convection with 
Robin thermal BCs of fixed Biot number 77 at the fluid boundaries is treated analogously. 

The details of the calculations for these two cases differ at various points, when it is 
necessary to consider the contributions of the plates on the one hand, or of boundary 
terms on the other. At the level of the strong admissibility criterion for background fields 
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and of formulae for the bounds for piecewise linear backgrounds with a given i5, however, 
the problems with and without plates map onto one another when rj = a = d/X: 

As pointed out in 14.31 the strong admissibility criterion on a background f{z) for the 
full plate-fluid-plate system, Qf,R^[v,6] > for Qf,R^ defined in (|4.2ip . depends only 
on the restriction of f onto the fluid domain z S [0, 1]. Consequently, for a given Re it 
coincides with the strong admissibility criterion of Appendix IB . 2 1 on r(z) for convection 
with thermal BCs applied to the fluid boundaries, Qt,r^ [v, 0] > for the quadratic form 
Qr.fie from (jB 12p . This is because in both criteria we optimize over the same classes of 
fields V and 9, as no BCs on 6 are assumed. 

Of course the background fields with and without plates, f{z) and t(z), should satisfy 
their appropriate thermal BCs, (|4.2p - (|4.3p or (|2.17p . respectively. However, for plates 
with a given d and A, (j5.3p shows that the piecewise linear background fg, defined in 
(|5.ip and satisfying (|4.2p - (|4.3p . automatically also satisfies mixed BCs with Biot number 
ri = a = d/ X. Thus for a given (5, over the fiuid domain < ^ < 1, f5(2;) coincides with 
Ts{z) defined in (IB 23p (so that also 7 = 7, Ar = At). That is, for a given Re a piecewise 
linear background fs is strongly admissible in the sense of S I4.3l (this is a condition on S) if 
and only if its restriction ts = fs\[o^i] is strongly admissible in the sense of Appendix IB. 21 

Furthermore, comparing (|5.5p - (|5.6p with (jB 26p - (|B 27p . the corresponding bounds due 
to strongly admissible piecewise linear backgrounds at a given S agree: Bpv,i.cr{S,b) = 
Bp^i,r,{S, b) and Ppwi.o-Ci^, h) = 2^pwi,r,('5, b) for tr = 77. 

In this analysis we have thus systematically mapped the conservative bounding prob- 
lem with imperfectly conducting plates onto that with mixed BCs with the fixed Biot 
number 77 = d/X. In the following sections, we discuss bounds for convection for these 
two problems simultaneously, assuming rj — a — d/X. The results are presented mainly 
in the notation of Appendix [Bl using rj, Tg, ^ and Ar, recalling that for piecewise linear 
backgrounds with the same 5, we also have 7 = 7, Af = At, and Tal^o.i] = Tg. 



5.3. Cauchy-Schwarz estimates on the quadratic form 

Recall the strong admissibility criterion for the background field ts{z): QTs,ReW,0] > 0, 
or in Fourier space (by (I4.24p - (l4.25p ) Qk = Qk;ra,fl:Jwk, ^k] > for all k and for all 
sufficiently smooth (complex- valued) functions Wk{z), 9k.{z), where satisfies Wk = 
DzDk = at z = 0, 1 while no BCs are assumed for 6'k. (However, thermal BCs enter 
the strong admissibility condition through the BCs for rg, which fix the value of 7 for 
given S and 77.). For piecewise linear background fields Tg{z) of the form (|B 23p (or fg{z) 
as in (|5.ip ). this criterion reduces to a requirement that S is sufficiently small, for given 
Re = bR/{b-l). 

Elementary Cauchy-Schwarz and Young inequalities applied to the Fourier space quad- 
ratic form Qk allow us to derive explicit sufficient conditions on 5 so that Qk > for 
all k, and hence to estimate upper boimds on Nu. To do so, we need to control the only 
indefinite term in Qk, /J^ 2r^ Re[wk^t], by t he ot her terms. For completeness we review 
the necessary estimates from lOtero et al. (|2002l) : Since i&k and Di&k (and hence also 
w\i ^k) vanish at both boundaries, we have 



\M^^)et{z)\^ 



^[wM dc 



< / |u-kD^i^|rfC+ / K^wu\dC, (5.7) 



where for < z < ^, by the Fundamental Theorem of Calculus and the Cauchy-Schwarz 
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kk(^)| 
|Du>k(^)| = 



D^wk dC, 



< 



< 



^z(^j^ \Dw^{0\' ddj < V^llDziklljgi,, (5.8) 
V~z(^j^ \Q^w^{0\^ d^y < V^llD^u^klljg ij. (5.9) 



Substituting these estimates into (I5.7P and again applying the Cauchy-Schwarz in- 
equality, for < z < i we obtain 



\w^.{z)ei{z)\< CdC 



1/2 



|D^-.||[oi, 



iDflkl'dC 



1/2 



1/2- 



< 



fli Dwk 1 

[0,2! ai 



-ilD^kll^i +f|||D2u;k||? 1 +-||^k|P„i 



fc2" """"[0,5] ' 02" ^"[0,2J 

(5.10) 



where we have also applied Young's inequality pq < \{ajp'^ + /aj) for any aj > 0. 
Proceeding similarly, we obtain an analogous estimate for ^ < z <1. 

For the piecewise linear background ts{z), for which r^- = — 7 < for < z < ^ and 
1 — (5 < z < 1, and t'^ = Q otherwise, applying these estimates we have 



< 





.5 


il 


) 






4^/2 


ai 



l^kfi'kl dz 



ai||Du)k||rr 



[oaj 



l-<5 
1 

ai 



|d4||[o,ij I p 



^^2 ||„2 - ||2 



a2 



^kllfoa] 



and thus 

r-l 



2r;Re[uik6'k] dz ^ j r'g (wJl + wl§k^ dz 



> 



ai||D^k|P + -||D4|P + ||||D2u,k|P 
2V2 L oi A;^ 



a2 



(5.11) 



where norms are taken over the entire interval [0, 1] unless otherwise indicated. Substi- 
tuting this estimate on the indefinite term into Qk given by (j4.25p . we find 



- ' i?c 2V2 

752 



1 jS^a2\ 1 



1 



-llD'^klP + ^fc^ll^ikir 



1 



2V2 



02 



fc'ii^kir 



2V2ai 



|Df?k| 



In the absence of any additional a priori information, fo r instance on the deca y rate 
of the Fourier coefficients (compare IConstantin fc Doerinel ( 1996[ ): Kerswelll ( 200lh ). our 
remaining estimates are necessarily fc- independent; we ensure the positivity of Qk by re- 
quiring all coefficients to be nonnegative. We choose 01=02 = 7J2/2\/2; then, dropping 
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manifestly nonnegative terms, 



We can thus guarantee that Qk > independent of k (and hence that ts is strongly 
admissible) if we choose 7^5^/8 < 1/Rc- For given thermal BCs, 7 = -f{6) is specified as 
a function of S; so this is a constraint on S to have Qk > 0, that is, for ts{z) to be an 
admissible background. Defining Sc by 

^^^^^ ^ 4is:TW '"wT' ^'-''^ 

we obtain the best bound in this approach by choosing S — 6c', the piecewise linear profile 
Ts (or fs) is strongly admissible for any S < 6c- 



6. Explicit asymptotic bounds for convection with thin, highly 
conductive plates or mixed thermal boundary conditions 

Using piecewise linear background profiles and the estimates in § 15. 3[ we may now 
derive explicit analytical bounds on the growth of the Nusselt number Nu with the 
control parameter i?, and hence with the Rayleigh number Ra. 

The results are described below mainly in terms of the mathematical idealization of 
mixed (Robin) thermal BCs with Biot number rj, showing that one may interpolate 
between the fixed temperature (Dirichlet) and fixed flux (Neumann) limits in a unified 
formulation. However, as discussed in § 15.21 all results apply also to the more physical 
problem of a convection in a fluid bounded by imperfectly conducting plates of flnite, 
nonzero (scaled) thickness d and conductivity A, when rj — a G (0,oo). We shall remark 
on possible interpretations of our results for convection with plates when appropriate. 

In this Section we summarize the main asymptotic bounds; more details, including 
improved values of the prefactors obtained by numerical solution of the optimization 
proble m for piecewise linear background proflles, are given elsewhere (jWittenberg fc 



The asymptotic analytical and the numerical bounds obtained using piecewise 



linear background functions differ only in their prefactors; the scaling with respect to R 
and rj (or a) is the same in each case. 

We begin by reviewing the results for Dirichlet (77 = 0) and Neumann (ry = 00) BCs, 
since in the general case, depending on the relative sizes of 6 and rj, the scaling behaviour 
agrees with one or the other of these extremes. In fact we shall see that for any f] > 0, 
the i? — )• 00 asymptotic scaling behaviour is as in the flxed flux case. 

6.1. Fixed temperature boundary conditions 

In the case of Dirichlet BCs (77 = or cr = 0), we have AT = Ar = 1, R = Ra, and 
(|B 25|) implies 7 = 1/26. Thus the sufficient condition (|5.13p on 6 simplifies to 6 < 6c 
where 

a2 32 6-1 



One c an show that the optimal choice of b in this formulation is 60 = 3/2 (see lWittenberg fc Gao 



( 2010( )). for which i?c = 3i?, and hence 6 < 6c = 4-^/2/3 R ^1"^ is sufficient to obtain a 



rigorous bound. Since for this b = bo, (|B 261) becomes 



Nu = p< 6pwi.o(<5, f'o) = l-bo + ^ = -^ + ^, (6.2) 
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for any 5 < (5c, the best rigorous analytical bound on the Nusselt number using this 
approach, valid for all R sufficiently large that 5c < 1/2, is 

where we used the fact that for fixed temperature BCs, the control parameter R is the 
usual Rayleigh number Ra. 



6.2. Fixed flux boundary conditions 

In the opposite extreme, for Neumann BCs (jj = oo), we have /? = 7 = 1, so Ar = 26 
from (|B 24[) . and we bound AT from below using (jB 27p . Since b > 1, in order for the 
lower bound I'pwi,oo(i5, ^) = 1 — b + 2 d b on AT to remain positive as i? 00 ((5 — >■ 0), 
we need b — 1 = 0{6). Thus following lOtero all (120021 ) we choose 6 = 1 + cS and let c 
take its optimal value Coo = 1/2, so that the bound on AT becomes 

Nu-^ = AT > PpwLoo(<5, 1 + CooS) - 1 + (1 + 5/2){25 -1)^^ + 5^ ^ h. (6.4) 

The condition on S is as usual 6 < ^c, where with 7=1 and b = 1 + (5/2, the equation 
(|5.13p satisfied by 6c takes the form 



i^R-^^^ 
l + d/2 R 



(6.5) 



for large R, for which (5 -> 0; and hence Sc ^ 4^/'^i? Thus we have (using p.3p ) 



AT > Ppwi,oo(^c, 1 + <5c/2) - ^(5e - T^-R"'/', 



and so 



i?a = i?AT > i?I?pwLoo((5c, 1 + (5c/2) 
21/3 



3 _3 

2 21/ 
3 



TVm < - — i?i/3 < 
3 



Ra 



21/3 



1/2 



i?2/3, 



(6.6) 



as m 



Otero et all (2002). Note the scaling Nu < CiR^^^ in terms of the control parameter 



1 /9 

R, which translates to the usual scaling Nu < C2Ra ' . 



6.3. Mixed thermal boundary conditions 

For general mixed (Robin) thermal BCs with fixed Biot number 77 (or equivalently, for 
plates with nonzero, finite a = d/X), we need to estimate both AT and f3, using (jB 26p 
and (jB 271) . where 7 and Ar are given in terms of rj and S by (IB 25p . The sufficient 
condition S < Sc for ts to be (strongly) admissible, derived via the Cauchy-Schwarz 
estimates of 15.31 is that Sc satisfies (j5.13p . which (substituting for 7 from (|B 25P ) here 



takes the form 



(l + 277)%4 8 5-1 



4((5 + ?7)2 Rc b ^ ^ 

We shall see that in this general case with < r/ < 00, the scaling of the bounds depends 
on the relative sizes of i5 and 77, behaving either as in the fixed temperature limit (for 
5 > 77) or the fixed flux limit (for S < rj); but that for any 77 > 0, the asymptotic scaling 
properties as i? — >■ cxi are as for fixed flux boundary conditions: 
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The fixed temperature problem rj ~ as a singular limit: 

Recall that for Dirichlet thermal boundary conditions ri — 0, we have AT = At = 1, 
so that we obtain an upper bound on Nu for any b > (there is no concern that the 
lower bound 2?pwi,o on AT may become negative), and we can choose 6 — 1 = 0(1) for 
all 6. In this case = 0, though, 7 = 1/2S is not bounded above as i? — > 00 (5 — > 0), and 
hence neither is /?; the growth in the (upper bound for the) Nusselt number Nu = /3 in 
the fixed temperature case with increasing control parameter R = Ra is due to that of 
the (non-dimensional) boundary heat flux. 

The situation is quite different for any nonzero Biot number 77: since Q < 5 < 1/2, we 
have < (1 — 25)/2{5 + ry) = 7 — 1 < l/2r/, so that for each 77 > 0, now 7 is bounded 
above as (5 — 0. From (IB 26|) (and choosing 6 < 2) it follows that for all 77 > 0, we have 
the rigorous (though presumably weak) upper bound on the boundary heat flux 

/? < 1 + A < 1 + i ; (6.8) 

that is, /? saturates at a finite value as i? — !■ 00. On the other hand, the (non-dimensional) 
averaged temperature drop across the fluid AT is not bounded below away from zero: 
AT — > 00 Hence asymptotically for large R, the growth in the Nusselt number Nu = 
/3/AT (and in the corresponding bound) is due to the decrease in AT, rather than due 
to growth of /3. That is, for any 77 > the (asymptotic) behaviour and scaling is as in 
the fixed flux case; the flxed temperature problem is a singular limit of the bounding 
problem. (A simila r observation was made in the context of horizontal convection by 



Siggers et al\ (|2Q04I ).) 



Scaling for poorly conducting boundaries: 

The nature of the Nu-R scaling depends on whether 5 > 77 or (5 < 77, and hence on the 
value of 77. For sufficiently large Biot number (largely insulating boundary) 77 > 1/2, we 
always have S < rj. Since for such 77, 7 is approximately constant (1 < 7 < 1 + I/277 < 2; 
compare 7 = 1 for 77 = cxd), we see from ()6.7|) that a sufficient admissibility condition 
for Ts is 6 < Sc = 0(Rc ^^*), as in the flxed flux case. We choose b = \ + c5 < 3/2 for 
some c < 1, so /3 < 1 + 6(7 — 1) < 5/2, and there is no transition in scaling regimes; as 
in the fixed flux case, for all sufficiently large rj the growth in Nu with increasing R is 
due to the decrease in AT. This conclusion equivalently holds for thick and/or poorly 
conducting plates for which a = d/X>0{\). 

Scaling regimes for highly conductive plates: 

For relatively small Biot number (largely conducting boundary) rj < 1/2, on the other 
hand, it is possible to have S > rj for low enough thermal driving, thereby allowing 
for different scaling behaviours. In particular, we consider the case of small Biot number 
(77 <^ 1, near the fixed temperature limit). This is relevant (by the correspondence a = rj) 
to convection in a fluid bounded by conductive plates in the physically relevant limit of 
< (T ~ d/X <C 1; some implications for that situation are discussed in § 16.41 

As the control parameter R increases, one observes a transition between two distinct 
scaling regimes: 

Low Rayleigh numbers: the "fixed temperature" limit: 

f Recalling the nondimensionalization, observe that this does not imply that the dimensional 
averaged boundary heat flux is uniformly bounded above, or that the dimensional averaged 
temperature drop AT* decays to zero as R ^ 00. 
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For sufBcicntly small R, we have 5 ^ rj; in this limit, we find 7 ~ 1/26 and At « 10 
and the sufficiency condition (|6.7p is 5 < (5c = 0{Rc )■ Since Ar is bounded below 
away from zero, so is the lower bound I?p„i ,,((5, 6) = 1 + 6(Ar — 1) > 1 — 6/2 on AT 
for any fixed b < 2. Thus we may obtain a bound on Nu in this regime by choosing 
any 5 € (1,2), and by comparison with the fixed temperature problem, it is sufficient 
to choose 5—1 = 0(1), so that the "effective control parameter" — bR/{b ~ 1) is 
proportional to R, and Sc = 0{R~^^^). 

It follows that the bounds on /3 and AT scale as /3 < Bp^\,n{5,h) = 1 + 6(7 — 1) = 
0{5-^) = OiR^'"^) and AT > Ppwi,,, = 0(1); hence the growth in the Nusselt number 
is due to the growth in the dimensionless averaged boundary heat ffiix, Nu ~ /3/AT < 
A/'pwL,, = 0{5-^) = 0(i?^/2). Furthermore, we have Ra = i?AT « i?, so that the control 
parameter approximately coincides with the usual Rayleigh number in this case, and we 
have 5 ^ 0{Ra-'^/^), and Nu < CiRa}/"^ for some 77-independent constant Ci. Hence 
when 77 ^ 1, for sufficiently small i?, everything scales as in the fixed temperature case. 

Transition: 

As the control parameter R increases, 5 shrinks, eventually decreasing below the Biot 
number 77; the system then enters another scaling regime, in which the above estimates 
no longer apply. The transition at (5 w 77 occurs (based on the low-i? "fixed temperature" 
scaling, which gives 5 ^ 0{R-^/^) ^ 0{Ra-^'^) in our formalism) for 

i? « i?t = 0(77-2), (6.9) 

that is. Rat = 0{ri''^). 

High Rayleigh numbers: "fixed flux" scaling: 

Once the "boundary layer thickness" 5 has decreased below 77 > for increasing R> Rt, 
we enter another regime (which does not exist in the fixed temperature case 7; = 0), with 
changes in the scaling behaviour of the bounds, and especially in the relative contributions 
of /3 and AT to the Nusselt number. 

In this regime, as i? — >■ 00 (that is, Ra — > 00) for fixed 77 the growth in 7 saturates, 
while At = 0{6) decreases. Asymptotically ior S <^ ij < 1/2, we have 7 ^ (1 + 2r])/2r] = 
JmsLxiv)) while At ^ ^(1 + 2ri)/r] ^ 1, and for each fixed 7; > the behaviour is now as 
if we had Neumann thermal BCsU 

More generally, for < r/ < 1/2 and decreasing S < ri, we have 7 = 0(7;"^) and 
At = 0{S/ri). Consequently, in order for the lower bound I?pwi,r; = I — b + bAr on AT 
from (jB 27p to remain positive as (5 — 0, the so far arbitrary parameter b > 1 must 
be chosen as 6 = 1 + 0{S/r)). We then find that /S < Spwi,^ = 1 + 6(7 - 1) = 0{r]-'^) 
saturates, while AT > 0{5/t]); hence the growth in Nu is now due to the decay in AT. 
In this regime the scaling behaviours are Ra > 0{SR/ri), Rc = 0{riR/S) and 

S = O (77i/2i?-i/4) ^ O (ri^'^R-^/^) = O (Ra-^/^) ; (6.10) 

more precise asymptotic statements are given below, while implications for convection 
with plates are in § 16.41 

Asymptotic scaling of hounds for < 77 < 00; 

Having outlined the behaviour in the different regimes, we now derive the scaling of 
the bound on Nu in the limit of large driving, R — >■ cxi, so that 5-^1 and S <^ rj; 

X Proceeding more carefully, for 5 > r;, we have 1/45 < 7 = (1 + 2i])/2{5 + 77) < 1/5 and 
1/2 < Ar = 5(1 + 2?7)/(5 + 7?) < 1. 

t More precisely, for 5 < 7?, we have 7max/2 = (1 + '2ri)/4r] < 7 < (1 + 2r])/2r] = 7max, and 

(57max = (5(1 + 27?)/277 < At < (5(1 + 2ti)/7] = 2(57niax. 
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Witte nberg fc Gaol (|20in[ ) for a comparison with the optimal solution for piecewise 



linear backgrounds ts{z). (As usual all these results carry over directly to convection 
with plates for a — rj.) 

In the light of the previous discussion, to ensure a positive lower bound on AT for 
<5 <C we must take h = 1 + c5, where the optimal value of c turns out to be 

l + 2?7 



4?7 ■ 

Using this optimal choice of 6, the lower bound (jB 27|) on AT becomes 

/ (I + 277) 
5 ^rj 

^ 3 (5(l + 2?7) 

' 4 



(6.11) 



AT > Ppwi,„((5, 1 + Cr^S) ^ -c,j5 + (1 + c„(5)- 

5(1 + 277) 3 + 2(5 



■n 



(6.12) 



while similarly, the upper bound (jB 261) is 



/3 < 'Bpwi,,,((5, 1 + c^(5) = -c,,(5 + (1 + c^(5) 



l + 2?7 

2{5 + r^) 



l + 2?y 
2(<5 + v) 



1 + 277 
277 ' 



so that an upper bound on the Nusselt number for admissible i5 ^ 77 is 



Nu 



P <: Kf 1 ^ 1 4 + 5(l-25)/77 



2 

36' 



(6.13) 



(6.14) 



3 + 2,5 

compare (|6.2p and (|6.4p . 

We remark that the width (5bl of the thermal boundary l ayer is often related to the 
Nusselt number via Sbl = (2Nu)~^ (jNiemela &: Sreenivasan (l20066h 'l: our high-i? result 
for the piecewise linear background, 6 ~ {3Nu/2)~^ for 77 > (or (5 ~ (4A''m/3)~^ for 
77 = 0), may be interpreted as a systematic statement of such a boundary layer model. 

Returning to the computation of asymptotic bounds, we note from (|6.1ip that for 
77 > 1/2, = 1/2 + 1/477 < 1, while for rj < 1/2, c^6 = (1 + 2ri)S/Ari < 5/2r], so 



that whenever S <C min(77, 1) we have c^(5 <C 1; consequently b — 1 + CriS 
Rc = bR/{b — 1) ~' R/cjjS. In this case the condition (|6.7p is thus 



5^ = 32 



32- 



1 + 277 



SR- 



-6R-\ 



or Sr 



(l + 27;)2 (1 + 277)2 4?7 "1 + 27/ 

27^1/3(1 277) ^^/"^i?^^/"^. Substituting into the above bounds, we have 

3"^ 

3 1 
4~ 



iVT/ < 7Vpwl,,,(<5c, 1 + CrjSc) 

Ra = RAT > i?2?pwK^(5c, 1 + c^S,) 



1/^1 + 277 
3 



»7 



2r; 



ScR' 



V 



R^'\ 



1 + 277 



1 and 



(6.15) 



(6.16) 



77 



2/3 



(6.17) 



so that we obtain a bound on the asymptotic scaling as i? — >■ 00 of the Nusselt number 
with the Rayleigh number whenever r/ > 0: 

1/3 / ^ \ 1/3 



Nu < 



1 / 1 + 277 

7? 



77 



1 + 277 



Ra^'^ 



(6.18) 



independent of the Biot number. Observe in particular, by comparison with ()6.6p . that 
the prefactor \j2j21 is the same as for the fixed flux problem. 
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6.4. Heat transport in thin, highly conducting plates 

It is instructive to view the above scaling results in the experimentally realistic context of 
conductive plates with small, but nonzero, thickness hg and/or large, but finite, conduc- 
tivity As, relative to the properties of the fluid. In dimensionless terms, this corresponds 
to fixed small, nonzero cr, since we have < d = hg/h <C 1, 1 <C A = Xg/^f < oo, and 
thus < cr = d/X < 1. 

Observe that in this case we have 1 + 2cr w 1, so that (from (I2.22p and (j2.26p ) Q « 
A* = T* — T*, that is, temperatures are scaled with respect to the applied temperature 
difference across the entire system. This allows us to interpret the control parameter 



agh^ agh^ A* agh^ 



R=^e= ' ^^{t*-t:) (6.19) 

VfKf UfKfl + Za VfKf 

as the Rayleigh number measured in terms of the imposed temperature drop across the 
full plate-fluid-plate system, instead of the temperature difference across the fluid only. 

Letting J be a measure of the thermal boundary layer width, for sufficiently small R we 
have (5 3> cr; that is, the thermal boundary layer thickness, approximated by /i(5, is much 
greater than the plate thickness scaled by the conductivity ratio, given by hgXf /Xg. In 
this limit we have AT « 1, which implies that essentially the entire temperature drop 
across the system occurs across the fluid (and that i?a « i?). That is, for cr ^ 1 and 
sufficiently small driving, the thermal behaviour of the fluid is essentially unaffected by 
the presence and finite conductivity of the plates, and the commonly used approximation, 
that the fluid boundaries are held at fixed temperature, is appropriate. 

As R increases, 5 decreases, until eventually 5 ~ a; this occurs for R k, R^, where in 
our analysis the transition value i?t — 0{a^^) (see (|6.9p ). Near this transition, Nu < 
0(i?V2) ^ e)(f7-i) = 0{X/d) (with 0(1) constant prefactors), so that we can interpret 
the scaling transition as occurring when the effective conductivity of the fluid, measured 
by the Nusselt number, becomes comparable to the plate-fluid conductivity ratio, scaled 
by the plate-fluid thickness ratio; this is in accord with our intuition. 

Once the control parameter R increases beyond i?t, the high- i?a asymptotic regime is 
entered, in which the scaling behaviours of the bounds differ from those in the low-i?a 
case. In particular, in the i? — >■ oo limit, when (5 <C cr ^ 1, using (|5.5p - (|5.6p and the 
^-scaling (|6.10l) we may estimate the bounds on AT and /3 in this analysis to be 

AT > Pp^i,, = O (a-^S) - O (a-2/3^-i/3A ^ ^ < ^^^^^ ^ ^ (^-i) ^ (g 20) 



and the usual Rayleigh number is related to R via Ra = i?AT > O (a 2/3^2/3~j jg 
apparent that for fixed nonzero cr, for suf ficiently large Ra al l the intermediate variables 



scale as in the fixed flux case discussed in lOtero et al\ (|2002l ). as expected 



In this scaling regime, the dimensionless averaged heat flux /3 through the fluid bound- 
aries saturates, while an appreciable portion of the temperature drop across the system 
now occurs across the plates, whose finite thickness and conductivity become significant 
for R> Rt. Consequently Nu increases no longer via growth in j3, but due to the decrease 
in the averaged temperature drop across the fluid, as a fraction of the overall applied 
temperature drop, according to AT > O (j^R/ Rt)~^^^^ ■ Finally, we flnd the high-i? 

asymptotic scaling of the bound on the Nusselt number Nu — /3/AT < 0{d^^) for con- 
vection in the presence of conductive plates with < ct ^ 1 (in this formalism, using 
a family of piecewise linear backgrounds and conservative Cauchy-Schwarz estimates): 
Nu < AApwi,. = O (a-i/2i?e'/') = O (a-i/3i?i/3) = o (iZa^/^). 

In summary, for Rayleigh-Benard convection in a fluid bounded by thin, highly but 



30 



Ralf Wittenberg 



not perfectly conducting plates, the main analytical results are: there exist cr-independent 
constants C2 and C3 so that as i? — >• 00, 

Nu < C2i^a^/^ (6.21) 
iVw < CaCT^^/^i?^/^ (6.22) 

The scaling in (|6.2ip is the same as obtained elsewhere for finite Prandtl number 

Rayleigh-Benard convection; in this formalism, the presence of conductive plates does 

not appear to alter the asymptotic scaling dependence of the Nusselt number on the 

usual Rayleigh number Ra. 

In contrast, consider the result (16.221) : while this (possibly non-optimal) bound on Nu 
1/2 

scales as Ra ' in terms of the Rayleigh number measuring the averaged temperature 
drop across the fluid, for sufficiently large imposed temperature gradient we find that Nu 
scales as R^^^ in terms of the Rayleigh number measured across the entire system; albeit 
with a prefactor that grows for small a as a~^^^. 

In an experiment with sufficiently small fixed dimensionless plate thickness d and/or 
large conductivity ratio A, so cr = d/X <C 1, it might seem plausible to ignore the plates 
and evaluate the Rayleigh number assuming that the fixed temperature difference is 
imposed at the boundaries of the fluid. We have shown directly from the governing PDEs 
that in terms of this "Rayleigh number" R across the full system, for sufficiently strong 
heating the scaling exponent p in a relationship Nu ~ CR^ could be no greater than 1/3. 
We should emphasize though that this "1/3 scaling" in our bounds is only relevant for 
large R (or Ra) — beyond a transition value i?t which scales, in our estimates, as cr^^, and 
may thus for small a be inaccessible to experiments or direct numerical simulations — 
and is presumably unrelated to the exponents p ^ 1/3 seen in experiments or simulations. 

7. Conclusions 

For finite Prandtl number Rayleigh-Benard convection, we have formulated the energy 
identities and bounding problem in the case of mixed thermal BCs with fixed Biot number 
T] applied at the upper and lower boundaries of the fiuid, and demonstrated that the fixed 
temperature and fixed flux extremes may indeed be treated as special cases of a more 
general model, for which one can obtain rigorous analytical and asymptotic bounds on 
convective heat transport. 

It has also come out of this formalism that, at least at the level of our conservative 
upper bounds with piecewise linear backgrounds, the case of convection with plates may 
be systematically mapped onto that with flnite Biot number, via rj = a = d/X. 

While the scaling of these analytical bounds on the Nu-Ra relationship remains well 
above that observed experimentally or in direct numerical simulations, some of the qual- 
itative conclusions may be instructive. Of particular interest is that while for each fixed 
control parameter R the bounds depend smoothly on rj for < ?? < 00, the asymptotic 
i? — 7- 00 behaviour of the bounding problem for any nonzero Biot number is as for the 
T] = 00 fixed flux problem. Indeed, we have proved that unlike in the fixed temperature 
case 77 = 0, for each ry > the averaged dimensionless boundary heat flux /3 is bounded 
above uniformly in i?, /3 — 1 < 77"^, so that the asymptotic growth in Nu is necessarily 
due to decay of AT. That is, the limits 77 — and i? — 00 do not commute: the much- 
studied fixed temperature case is a singular limit of the general bounding problem. From 
the point of view of understanding a realistic convection situation in the limit of large 
R, it appears that the mathematical structure of the insulating-plates fixed flux problem 
is more relevant. 
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Furthermore, our analysis reveals two distinct scaling behaviours for sufficiently small 
nonzero cr (or 77): In the "fixed temperature scaling regime" for small Rayleigh number Ra, 
the growth in the (bounds on the) convective heat transport measured by Nu is largely 
due to the increase in the averaged boundary heat flux /3. However, for strong driving, 
eventually a "fixed fiux scaling regime" is reached in which the effective conductivity 
of the fluid due to convective transport exceeds the plate conductivity, and the plates 
effectively act as insulators; /? saturates and further increases in Nu are due to decreases 
in the averaged temperature drop AT. The transition between these regimes occurs when 
the "thermal boundary layer width" 5 is comparable to cr. 

It would be of interest to determine whether this qualitative transition at i? « _Rt from 
effectively conducting to effectively insulating boundaries is in fact reflected in the physics 
of convective turbulence in the fluid, and thus observable in experiments with small 
a — d/\^ or in direct numerical simulations with fixed Biot number 77 ^ 1. It is possible 
that it may not be: the re cent direct numerical s imula tions comparing fixed temperature 
and fixed flux BCs, due to Johnston fc Doerind ( 20091 ) in two dimensions with horizontal 



periodicity, and to IVerzicco fc Sreenivasa d (|2008|) and IStevens et all \2Q1&] in three- 



dimensional cylindrical geometry, suggest that the heat transport in large- i?a turbulent 
convection appears to be insensitive to thermal boundary conditions. 

In this context we observe that the prefactor in our asymptotic analytical bound Nu < 
CRa^/"^ increases from Co = 3\/6/32 « 0.230 to = Coo = ~ 0.272 for 

77 > 0; that is, within the framework of our upper bounding calculations with piecewise 
linear background it appears that the estimates on the heat transport increase when the 
boundaries are not perfectly conducting. It remains to determine whether this increase is 
an artifact of the choice of background t(z) or of the background flow bounding approach 
in general. 
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Appendix A. Comments on the formulation and notation 

In the following we introduce and clarify some notation used in our calculations. 

Averages: 

Following Otero et al\ ( 20021 ). for functions h{x, y, z) and g{t) we define the horizontal 
and time averages, h{z) and {g) respectively, by 



and 



(5) = limsupi / git)dt, (A 2) 

r— )-oo ^ Jq 



where A = L^Ly is the nondimensional area of the plates. 
Global definitions in presence of plates: 

The problem of Rayleigh-Benard convection with bounding plates may be formulated 
using separate fields in the fiuid and the lower and upper plates, with appropriate condi- 
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tions at the interfaces between the different domains. However, as discussed before ()2.19p . 
for simphcity of notation it is more convenient to treat the space- and time-dependent 
fields as being defined across the entire plate-fluid-plate system, for —d < z < 1 + d. 
(Recall that all x- and y-dependent quantities are L^, Lj^-periodic in the horizontal di- 
rections.) Thus we consider a single dimensionless temperature field T, which coincides 
with the temperature in the lower plate Tp i on z G [— d, 0), with the fluid temperature 
Tf on z Cz (0, f ) and with the upper plate temperature Tp „ on z S (f , 1 -I- c?]; it is a 
continuous function with discontinuous vertical derivative at z = and 1, which sat- 
isfies the conditions (|2.30l) - (|2.3f p at the fluid-plate interfaces at z = and 1, and the 
boundary conditions (j2.32l) . (Equivalently, defining a piecewise constant global thermal 
conductivity function A — A(z) which takes the values f in the fluid and A in the plates, 
(j2.30p - (|2.3f p can be interpreted as continuity conditions on both T and the weighted 
derivative A dT/dz.) Similarly, we may extend the definition of the velocity field by u = 
in the plates —d < z < and 1 < z < 1 -I- d; then (by the no-slip BCs ul^^o.i — 0) the 
velocity field is similarly continuous across the entire system, with discontinuous vertical 
derivative in the horizontal velocity components (by incompressibility, w^l^^o.i = 0). 
Similar considerations apply to the fiuctuating quantities w — \i^9 = T— f defined in 

Limits and boundary terms: 

For convection in the presence of plates, since the temperature field T is piecewise 
defined, care should be taken in evaluating and related fields which are discontinu- 
ous at the interfaces z — Q and 1, for instance when evaluating boundary terms upon 
integrating over the fiuid or plates. To simplify the description, before (|2.19p we in- 
troduced notation for limits, writing, for instance, T\z=i+ = lim^^i+ T = Tp,u\z=i, or 
[dT /dz)\z=o+ = {dTf /dz)\z=o. Similarly, we write (•)|^Io-i- *o indicate that boundary 
values are approached from within the fluid: specifically, for any function /(z), we have 

(/)l':o+^/l-i--/l-o+^ lim /- hm /= T/.dz; (A3) 

z — ^1 — z — ^U~r iJ 

and similarly for (/)|°I_^ = fz dz and if)\lti+ = Ii^"^ fz dz. 
Integrals: 

For a function h{x, y, z), we define volume integrals of h over the fiuid, lower plate and 
upper plate by h, h and h, in the expected way: Over the full fluid layer, we have 

h = A h{z) dz — I 11 h{x,y, z) dxdy dz ; (A 4) 

// Jo ^0 JJa 

while over the lower and upper plates, respectively, 

h^A h{z)dz, = h{z)dz. (A 5) 

II J-d Ju Jl 

The usual norm is deflned over the fluid layer by 

1.2 



/ 



hWj^ I ^ I II h^{x,y,z)dxdydz. (A 6) 



For the full convection problem with plates, it turns out that the energy identities are 
best formulated in terms of a conductivity-weighted integral across the plate-fluid-plate 
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system; we define 

{h} = X h+ h + X h= / / h{x, y, z) dxdy A(z) dz = A I h{z) A(z) dz, 

Jl Jf Ju J-d J J A J ~d 

(A7) 

which may be interpreted as an integral with a weighted measure d(^ = X{z) dz. 

With this notation, the divergence theorem applied to vector fields H over the fluid 
gives /y. V • H = v4 H • Gz I ^^g_|_ , using horizontal periodicity, and similarly for integrals 
over the plates. 

Notation for quantities defined in presence of plates: 

A further notational convention we introduceU is the use of a bar ~ to denote quantities 
relevant to the problem with plates, such as the background field f(z) (defined on [—d, 1 + 
d]) with associated 7, Af and fa, and the bounds Ba, and (or, for piecewise linear 
backgrounds, Spwi.o-, 2^pwi,<T and A/'pwi,^)- This convention is chosen to distinguish them 
from the corresponding quantities (such as the background t(z) defined only on [0, 1]) 
for the convection problem without plates treated in Appendix |BJ in which thermal BCs 
are applied directly at the boundaries of the fiuid. 



Appendix B. Derivation of bounding principle for mixed (fixed Biot 
number) thermal boundary conditions 

In §§ I3.2H5.1I of the main body of this manuscript, we have chosen to concentrate on 
convection in a fiuid bounded by conducting plates of finite thickness and conductiv- 
ity with fixed temperatures applied to the outer boundaries of the plates, deriving the 
governing identities and bounding formalism for that case. In this appendix, in parallel 
with the presentation in the main text we obtain the analogous formulae for mixed ther- 
mal boundary conditions applied directly to the fiuid boundaries. In many cases we can 
reuse computations over the fiuid layer < z < 1 (in nondimensional variables), merely 
modifying the BCs on the temperature at z = and 1. 

We recall from § l2.2l the definition of the mixed (Robin) thermal BCs; specifically, from 
(|2.17p we have, in dimensionless form, that for a fixed 77 (we assume here < r/ < 00) 
the temperature field T satisfies T — ryT^ = 1 + r/ on z = 0, and T + rjT^ = —rj on z = 1. 

B.l. Governing identities 

Relation betiueen /3 and AT: 

As in § 13. 2| we can relate the averaged (nondimensional) boundary temperature drop 
AT and heat fiux /?: taking horizontal averages of ()2.17p . and subtracting the upper 
boundary condition from the lower, we find T\z=q — T|z=i — r/C^ z\z=o + T z\z=i) = 1 + ^rj. 
Taking time averages and using (j3.ip and (j3.7p . we find the fundamental relation 

AT + 2ry/3 = 1 + 2?/. (B 1) 

Hence for < 77 < cxd, an upper bound on j3 constitutes a lower bound on AT, and vice 
versa, and we only need to bound one of these quantities to obtain an upper bound on 
Nu = /3/AT = /3/[l + 2?7(1 - /?)] = AT^^ + (AT"! - l) /2r]. In the following, assuming 
77 < 00 we shall solve for AT using (|B ip to present the identities in a form, valid in the 
fixed temperature limit 7/ — > 0, that allows the derivation of an upper bound on /3. 

f There should be no confusion between this notation and the earlier use of h to denote the 
horizontal average of a function h{x,y,z) or h{x,y, z,t). 
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Global energy identities: 

We next express the general energy identities over the fluid layer, p.l7p and p.l9p . 
for fixed (finite) Biot number BCs. Substituting AT = 1 + 2r]{l - (3) into ((XTT)) . we find 

^(||Vu||^.) =/3-AT=(l + 2r;)(/3-l). (B 2) 

Similarly, we can evaluate the boundary term in p.l9|) by using the BCs (|2.17p to solve 
for r at z = 0, 1 in terms of T^; substituting into TT^I^^g and taking horizontal and 
time averages, we find that the global thermal energy identity p.lQp becomes 

\{\\VT\\]) = {TT:\\J = (1 + 2,7)/3 - r;(lf + Tf (B3) 
(note the additional quadratic boundary terms not present in p.21|) ). 

B.2. Background fields 

Following the "background flow" variational method for obtaining upper bounds, as in 
§ 14.11 we decompose the velocity and temperature fields via u(x, t) = v(x, t), T{x,t) = 
t{z) + 6'(x, t)^ The background field t{z) is assumed to inherit the BCs on the temper- 
ature; assuming that the upper and lower boundaries of the fluid have identical thermal 
properties, we require t'(0) = t'(1) (compare ([BTZ])), and define 

At = t(0) - t(1), 7 = -t'(0) = -r'(l). (B 4) 

In the case of fixed Biot number BCs, r(z) satisfies (|2.17|) . which implies that 

At + 2?77 = 1 + 277 (B 5) 

Consequently, the perturbation 9 satisfies the homogeneous Robin BCs 9 + rjn ■ V0 = 
at the interfaces, which in our geometry become 

9 ~ r]9, ^ at z = 0, 9 + rj 9, = at z ^ 1. (B6) 

Evolution equations and identities for fluctuating fields: 

Substituting the decomposition of u and T into (|2.6p - (|2.8p (for general thermal BCs 
on the fluid boundaries), we readily find that the fields v and 9 evolve in the fluid as in 
(|4.7p . (I4.8P and (|4.10p (with r instead of f), with homogeneous BCs. The evolution 
for the perturbed temperature 9, found by multiplying (j4.10p by 9 and integrating, is 

Comparing with (|4.12p . we observe that when general (non-Dirichlet) thermal BCs are 
imposed at the fluid boundaries, boundary terms remain and play a significant role. We 
also use the relations u = v,T = t + 9 between the physical fields and the fluctuations 
over which the optimization is performed, expressed in the form 

l||Vu||^. = i||Vv||^, (B8) 
II VT||^. = II V^ll^. + ^J^O,t' + r'2. (B 9) 

As in § |4?2l we now form the hnear combination 6 • [2 • (|B7)) + ([B9)) ] + (1 - 6) • [ (|B8|) ]. 
t See Appendix lAl concerning tlie use of a bar " to distinguisli between, say, t{z) and f{z). 
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Substituting (|B 2[) and (|B 3[) . this becomes 



where as in ()4.2ip we define the quadratic form 



J-bQrMA^,9], 
(Bll) 



—\Vvf + 2T'we + \V0{^ 

He 



(B12) 



with i?c = bR/{b-l) (sec gH])). 

Using T = T + 9, (jB 41) and the definitions (13.11) and p.7p . we may write (6't' 
— 7(6*1^ ^) = 7A0, and decompose the first term in (IB 11[) via 



- "1 ) 

2=0' 



iLo) 



(B13) 



Substituting into (jB lip , writing A0 = AT — At and rearranging terms, we obtain 

b {(3 At - J AT) + {1 - b){l3 - AT) = b (^^' r'^ dz - 7Ar^ - ^Q^Jv, 9]. (B 14) 

Here we have introduced a modified quadratic form with boundary terms added to (|B 12[) 
(these extra terms vanish in both the fixed temperature and fixed flux hmits) , 



Qr,Rd^,0]^Qr.nAv,9]-A 



1 1 

\z=0 



(B15) 



The identity (|B 14[) is independent of the thermal BCs apphed directly at the fluid 
boundaries. In speciahzing to mixed, fixed Biot number BCs (|2.17p . we solve for AT and 
At using (|B 1[) and (|B 5p (for 77 < 00) to compute 

/3Ar - 7Ar + 7Ar = (1 + 2j])P - 2t]j'^. 

Now substituting (IB 2p and (IB 161) , (|B 14^ becomes in terms of /3 



(B16) 



(l + 2,y)(/3-l) 



r'^dz - (1 + 277) + 2777^ 



(B17) 



where for fixed, finite Biot number, the boundary term in Q'^ from (jB ISP can be 
written, using (jB 6p . as 



,1 

1 2=0 



{0lU=o + 9^A 



(B18) 



Observe the similarities between (IB 171) and the analogous formula (j4.19p for convection 
in the presence of conducting plates; in this case, though, explicit contributions from the 
fluid boundaries replace the integrals over the plates. 

A bounding principle for mixed thermal boundary conditions: 

As in S I4.3j we can now use (jB 17P to deduce an approach to bounding (3, and hence the 
Nusselt number Nu. We begin by letting allowed fields be those sufficiently smooth scalar 
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fields 9 and divcrgencc-free vector fields v which satisfy the boundary conditions of the 
problem. A background profile t(z) satisfying the thermal BCs will then be admissible 
if, for a given i?o > 0, the quadratic form Q'^ ^ appearing in (jB 17p is nonnegative, 
Q'^ jj. [v, 9] > for all allowed fields v and 9 (and hence for all solutions of the evolution 
PDEs for the fiuctuating fields). 

We now formulate an upper bounding principle for the Nusselt number Nu in terms of 
the control parameter R (and hence in terms of Ra — RAT): For each i? > 0, if we can 
choose 6 > 1 and a corresponding admissible background field t(z), then from (|B 17p the 
averaged boundary temperature gradient (3 is bounded above (for 77 < oo) by 

P<l-b+ (^^ r'2 dz + 2,77') - B,[t; b], (B 19) 

while using ()B ip and (|B sp . a corresponding lower bound for the averaged temperature 
drop across the fluid AT is (for ?7 > 0) 

AT > 1 + &(2Ar- 1) -b '^^ ( [ t'^ dz + —AtA =VJT:b], (B 20) 

1 + 27] \Jo 2ri J 

where the above equations define the functionals S,j[r;6] and 2?^[r;6] (where I?,Jt; 6] + 
2ryi3^[r; b] = 1 + 2rj). It follows from p. lip that the Nusselt number is bounded above by 

Nu<Nn[T]b] =Bn[T:b]/Vn[T;b]. (B21) 

Observe that for the conduction solution t(z) — \ — z, we have Bj^[t; b] — V^jIt; b] = 1, 
so that whenever this is an admissible profile, the bound on the Nusselt number takes 
its minimum value of 1, as expected. 



Strong admissibility: 

As in § 14. 3[ for convection with fixed Biot number BCs it is convenient to strengthen 
the admissibility condition on r(z), here by ignoring the boundary conditions on 9 and 
removing the boundary terms from the quadratic form ^ : 

We say that a background field t(z) defined on [0, 1] is strongly admissible if the 
quadratic form Qr,R^ defined in (jB 12p is nonnegative, Qr.fleiv, 0] > for all sufficiently 
smooth, horizontally periodic v and 9, where V • v = and v = on z = 0, 1, but we 
impose no additional constraints on 9. 

That the strong admissibility condition (positivity of Qr,i?e) implies admissibility (pos- 
itivity oi Q'^ j^) for fixed Biot number boundaries follows in this case from the fact that 
by (jB ISp . the additional boundary term appearing in (jB 15p for 77 7^ 0, 00 is stabilizing: 

Q'rA ^] = Qr,R. 9]+A7](^ (dlU^o + = ) > Qr,R. [v, 9]. (B 22) 



B.3. Explicit conservative bounds using piecewise linear background profiles 

While the strong admissibility criterion is difficult to verify for general backgrounds t{z), 
by restricting consideration only to piecewise linear profiles we may obtain analytical 
conditions ensuring Qt.fi.A'^,S] > 0; such profiles shall in fact also allow us to map the 
convection problem with fixed Biot number BCs and that with plates onto one another. 
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Figure 4. The piecewise linear background profile ts{z), with r' 

and r' = in the bulk. 



-J in the boundary layer, 



Piecewise linear profiles for fixed Biot number conditions: 

We define a one-parameter family of piecewise linear background profiles ts{z) over 
the fluid domain z G [0, 1] as follows: 

r Ta-j{z-S), 0<Z<S, 
t(z) = Ts{z) ^ I Ta, d<Z<l-S, (B 23) 

[ r,-7(z-(l-^)), 1-S<z<l; 
see figureSl From (|B 23p we immediately compute r(0) = Ta +7(5, t(1) — Ta — 7(5, and so 

Ar = t(0) - r(l) = 2(57, / ^'^ = 2(57^ = 7AT, (B 24) 

Jo 

where it remains to choose the average = |(t(0) + t(1)) and boundary slope 7 = 
— r'(0) = — t'(1) of the background as functions of 6. Substituting At — 2Sj into (|B sp . 
we obtain the values of 7 and At (for given S and 77) for which ts{z) defined in (IB 23[) 
satisfies the mixed (Robin) thermal BCs: 

2[d + rjj o + r] 

Now using this 7 in the BC (|2.17p written as t(0) = 1 + 77 + 77t'(0) for rj < 00, we find 
that Ta = 1/2, completing the specification of the background ts{z). (In the fixed flux 
case ?7 = 00, in which the governing equations depend only on temperature gradients, t^ 
is arbitrary; for consistency we choose t^ = 1/2.) 

Substituting (|B24l) - (fB25l) into (|BT9| ~ (|B201) and simphfying, we now find that the 
conservative bounds on /3 and AT for fixed Biot number convection with a piecewise 
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linear (pwl) background profile ts take the concise forms 

/3<epwi,^(<5,6) 
Ar>2?pwi,„(<5, b) 

corresponding to an upper bound on the Nusselt number of Nu < A/'pwi, »; {S, b) = J\f,-i [rg ; b] = 

Spwl,r,(^, &)/2?pwhr,('5, b). 

Relation to convection with plates: 

The development in this Appendix of a bounding principle for convection with mixed 
thermal BCs of fixed Biot number 77, applied directly at the fluid boundaries, parallels the 
calculations in l3.2| - [5TT] for convection with imperfectly conducting bounding plates. 

As discussed in § I5.2[ it turns out that the strong admissibility criteria for convec- 
tion with and without plates coincide, so that the estimates of 15.31 for verifying strong 
admissibility for piecewise linear backgrounds also apply to the present mixed EC case. 
Furthermore, for convection with plates with a — d/X, the piecewise linear field fs{z) 
defined in (j5.ll) satisfies mixed thermal BCs with Biot number r] = a, and the corre- 
sponding bounds on Nu agree. Hence the further derivation of bounds for fixed Biot 
number BCs can proceed simultaneously with that for convection with bounding plates. 
Asymptotic bounds for convection with Biot number ry, and equivalently for plates with 
a — 7], are obtained in § [51 with a discussion of their scaling regimes for increasing R. 
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